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1.0 : SUMMARY - 

This report describes the effort of the LTV Aerospace Corporation 
Vbught Missiles and Space Company, Texas Division under Task 1.1.2 of NASA 
Contract NAS9-10R10. The purpose of this effort was to extend the. (vipability 
of on existing system of computer routines , the Radiator Design System (RDS), 
described in Reference 1, to produce ah Advanced Radiator Design System (ARDS). 

The ARDS consists of three separate routines which .are (l) thb steady 
state design routine, versioh-2 (sSDR-2), (2) the Transient Perfprmance Routine 
Version 2(TPR-2), Md(3) the Computer Incident Heat Data Routine (CIHR). 

The SS DR determines a weight optimized radiator (or set of radiator 
panels) and low load control technique which satisfy user input requirements 
for maximian and minimum heat rejection, manber and type of radiator panels, 
area limitations, structural siqjport weight penalties i fluid properties , 
radiator material properties, micrometeoroid protection requirements, pressure , 
drop (ahd associated" pumping power penalty ) j ininimumT allowable fluid temperature, 
and use of supplemental expendable heat, sink devices (such as water boilers 
or sublimators ) . It optimizes the radiator (s)-coritrol technique combination 
by variation of tube diameter* tube spacing, fin thickness and flow path 
arrarigement within limits specified hy the user. It then outputs a definition 
of all panel and control technique configurations (ihcluding regenerator weight , 
if applicable) within ten percent of k minimum' weight design. The latest 
version ( SSDR-2) can be used to design solar absorbers as well as radiators. 

In addition , the characterizations of two-dimensiohal fihs i. fluid heait transfer 
coefficients , friction factor -and Vei'ght- calc.ul’atiohs-.- are significantly im- • : ' 

'.■proved bver^the '■previous 'yersiori.’ ■^..■V " V: ■ ' ■ 

By use of the TPR, the user may select k SSDR design and quickly 
determine its transient perfomanc® a radiator subsystem consisting of ‘the 
radiator Euid heat load control valving, subsystem piping and/or a regenerative 
heat exchangers The results of the TPR cem be used'to Jndge thejadeqviacy of 
effective environments and effective heat loads used for the initial steady 
state design. : Tbe SSDR arid, TPR can be rim sequentially to iterate on a weight 
optimum. subsystem design which satisfies transient requirements.; The TPR 
includes automatic output; plotting, interruption and restart capability j pro- 
visions for use of input data from magnetic tapes which can be edited .to in- 
corporate design . changes , or changes for parametric per formabcs studies, The 
TPR is -facilitated by SUBANS and SUBRAD which generate nodal lumped parameter 
TPR input data from a user specification of byerall radiatbr panel dimensions, 
fluid flow path geomet^., nvukher of flow paths , and specification of zones 
for ils Crete .tiime yariarit incident heat yetlues . . ! SUEANS divides one dimensional ^ 
panels and SUBRAD divides radial panels. 

With these two subroutines, the user may ais 6 ' automatically generate 
nbdkl models •with varying degrees of. fineness to establish accuracy, require- 
ments consistent with preliminary ’design objectives and mihimization of 
■computer'' . time ■ us.ed.."''-*''' ■ 


The CIHR permits the user to specify the vehicle prhit and 
orientation and obtain the incident heat curves on punched cards in TPR 
format with the minimum number of points. No change was made to the CIHR 
under this contract. 
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2.0 INTRODUCTION 

The projected trends of future long diuration space vehicles 
indicate the need for minimum weight space radiator subsystems which can 
operate over heat load ranges that are much wider than the capabilities 
of current subsystem design. The upward trend in maximum heat rejection 
requirements of low temperature radiators required for environmental con- 
trol systems necessitates increases in radiator panel area to, the extent 
that deployed radiator panels in addition to multiple panels which are 
integral with the space vehicle structure may be required. However, during 
period of dormant or low activity operation, reduced space vehicle heat 
loads may cause complete radiator fluid stagnation or freezing due to the 
heat rejection capacity of large radiator areas. The limit of radiator low 
heat load operation (minimum of the heat load range) is established by the 
coolant fluid freezing point and by the heat load control method employed 
on the radiator panels (e.g., flow stoppage in a portion of the radiator 
tubes) and in the radiator subsystem (e.g. ’, fluid bypass, or regeneration) . 
The use of mechanical, shuttering devices to block the radiator view of 
the space environment is not considered practical for large area systems 
because of excessive weight and system complexity. The weight of large 
area integral and deployed panels must be minimized dxaring design since the 
radiators becbme a significant portion of the overall spacecraft weight. 

Experience with previous radiator systems has shown that transient 
performance is a key factor which must be considered during conceptual and 
preliminaiy spacecraft design phases to avoid defining excessive radiator 
area or over complicating low heat load radiator control techniques. The 
set of radiator design programs svanmarized herein were developed to provide 
the analyst with a capability to generate optimum weight radiator panels 
or sets of panels from practical design considerations, including transient 
performance. In addition, modifications to existing detailed transient 
radiator analysis programs to .improve capability and user convenience were 
implemented. 
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3.0 


STEADY STATE DESIGN ROUTINE 


The Steady State Design Routine Version 2 (SSDR-2) written to 
determine the optimum weight design of radiator or absorber panels is 
described in this section of the report. The SSDR is capable of steady state 
design of weight optimized radiator or solar absorber panels. consistent with 
given space vehicle operational characteristics and constraints. The panels 
designed with SSDR-2 can be easily analyzed with the Transient .Performance 
Routine (TPR) described in Section 4.0 to establish transient performance in 
a radiator subsystem (consisting of the radiator and heat load control 
valving and/or a regenerative heat exchanger). The analyst can use TPR results 
to iteratively rerun SSDR with modified design constraints in order to establish 
a weight optimized design which yields desired radiator subsystem transient 
response characteristics. 

3.1 GENERAL ROUTINE DESCRIPTION 

The SSDR-2 radiator or absorber design requirement^ are input into 
the routine by the user and include: the maximum heat load, minimum heat load, 

the inlet .temperature at maximum heat load, the required system exit tempera- 
ture (mix temperature), panel area availability, panel thermal properties, 
and the effective thermal environment at maximum and minimum heat load. The 
design requirements at maximiun- heat load are used to size the system for the 
various panel geometries considered. 

The SSDR-2 optimizes the system weight at maximum heat, load by 
variation of the tube diameter ,. tube spacing (number :of tubes), panel flow 
path connection (series and/or serpentine tube connections), and panel fin 
thickness. Each of these parameters are varied independently within limits : 
imposed by the program or input by the 'userl. The sy.stem configuration must 
then meet or exceed the minimum heat load design requirements- in order for it 
to be considered an optimum design. 

In addition tb the satisfying of the above design requirements, 
the radiator or absorber system optimized by the SSDR-2 must also meet or 
consider the following operational constraints:, the system panel pressure 
drop, must not exceed' the maximum allowed, the transport . fluid temperature at 
minimum heat load must not. be less than a specified input valued and the panel 
area used for heat rejection cannot exceed the maximum available. 

If the panel pressure drop exceeds the maximum allowed, the 
particular panel configuration analyzed is not considered for the optimum 
design. The maximum allowable pressure drop constraint is input data and may 
be deleted' by the user. In this case the SSDR..’2 utilizes a pumping power 
penalty (user input ) to account for the influence of pres surei’ drop in determining 
an optimum weight panel or set of panels. . 


1 System designs using previously developed, radiator or absor'ber paneir”corr- 
figurations can be analyzed by -user specification of one or all of these 
variable parameters.. The specified parameter overrides the instructions 
within the routine for its variation. 

2 'For solar, absor'bers, the fluid temperature at the minimum heat Toad must 

not be more than a specified value. 
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» The minimum allovable fluid temperature (user input) will force 
the SSDR-2 to reject all heat load control methods that result in radiator 
exit temperatures at minimum heat load operation, which are less than the 
minimum allowed. 

If the radiators cannot reject the maximum heat load with given 
radiator panel area, the routine assumes that an expendable heat rejection 
system is in the system to dissipate the excess heat. The excess heat dissi- 
pation with a boiler /sublimator is optional and may be deleted by the user. 

3.2 ANALYTICAL METHODS 

3.2.1 Groundrules and Assumptions 

Fifteen general SSDR-2 grotmdrules are tabulated in Table *3-1. 

The most significant of these is use of a pair of effective thermal environ- 
ments for steady state design at maximum and minimum heat load conditions. 

For many spacecraft missions the system heat load and thermal environment 
may be in transient modes for a majority of the mission duration. T^us the 
SSDR user must estimate a maximum effective absorbed heat (for example in the 
sunlit portion of a planet orbit) which will not result in oversized panel 
area. This is required because the panel will not reach equilibrium 
operating conditions at the peak environmenval heat absorption point in the 
orbit. Similarly in the shadowed portion of the orbit the radiator will not 
be at steady state conditions as it views the actual minimum environmental 
heat flux for a short period, leading to a need for an effective minimum 
absorbed heat input to SSDR-2 . If a .short duration minimxim environmental 
heat absorption is used, an overly complex low heat load control method may 
be specified by the SSDR-2. 

The transient nat\ire of the radiator heat loads (rejection require- 
ments) should also be considered when using the SSDR-2. Specification of 
equilibrivun heat loads which will not be realized doiring the mission can 
similarly lead to unnecessary complication of radiator design requirements. 

If a panel is to be designed to satisfy such transient mission ■. 
requirements, the user can estimate effective absorbed heats and heat loads 
initially, run SSDR-2, and use the routine design in TPR to evaluate the 
reasonableness of his estimate. The two programs can thus be run sequentially 
to iterate on an optimum radiator subsystem design. 

Multiple panel systems may be designed by the SSDR-2. 

These systems may consist of panels which are identical to each other or 
panels which may be of different geometries (width, tube diameter, fin 
thickness). The panels may be flowed connected in series or in parallel 
depending on user specification. Ttie effective thermal environment may be 
different for each panel. A detailed description of the constraints of 
multiple panel systems follows. 

Panels may be connected for parallel flow; with the following 

conditions: 
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TABLE 3-1 

GENERAL GROUNDRULES 

a. The SSDR-2 utilizes a. pair of effective absorbed heat fluxes to establish 
the radiator design under maximum and minimum environmental heat influx 
conditions. 

b. The radiator or absorber system consists of fin-tube panels and heat load 
control methods; 

c. The panels may be connected in series or in parallel as des- 
cribed in Section 3 .2. 2. 

d. The panels are fin-tube surfaces with simple rectangular fins 

_and round tubes^ _ _ _ _ . 

e. The panel edges are. adiabatic. 

f. Heat transfer resistance of the surface coatings is assumed negligible. 

g. Steady state operation of the system is assumed. 

h. The tube temperature is assiuned constant around the tube periphery. 

i. The tube edge of ,radiati.ng fin segnehts are assumed to have uniform 
longitudinal temperature,. 

J. The temperature across the fin thickness is assumed uniform. 

k. The fin midppint between flowing- tubes is adiabatic . 

l. The absorbed heat (or equivalent sink temperature) is uniform for a 
radiator panel. . 

m. Longitudinal conduction in the fin, tube, and fluid is neglected. 

n. The radiant interchange between the panel fins and tubes is 
assumed to be zero. 

o. The tubes and fins of all panels are of a uniform material. 

The material properties are evaluated at the panel base 'temperature. 
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a. Parallel flow systems are limited to two panels flowed 
in parallel. 

b. Paralle]. panels are of identical design. 

c. The absorbed heat may be different for each panel - 

d. The coolant flow is equally divided between the panels 
unless the sink temperature is greater than the inlet 
temperature for one panel. If the sink temperature is 
greater than the inlet temperature for one panel the flow 
bypasses that panel. 

Figure 3-l(a) shows a. typical radiator or absorber system with 
the panels connected for parallel . flow. 

tna,y tji with thb roilawing 

conditions: 

. . a. Panels connected in series are limited to 10 panels. 

b. Panels may be of two different types. Panels of the 

same type are identical to each other for all the variables 
being considered (e.g. fin thickness, tube spacing, etc.). 

c. All panels of the same type .must be sequentially arranged. 

d. The absorbed heat may be different for each panel. 

Figure 3-1 (b) shows a typical system with the panels connected 
for series flow. 

Fourteen heat load control- methods considered by the SSDR-2 to 
assure adequate system performance at minimum heat load ure described in 
Table 3-2. The control methods consist of bypass, regenerative, valve 
stagnation, and various combinations of all three of these. The control 
methods are in many instances combined with panel tube patterns to yield 
wide range heat load radiators. By inputting a zero for the minimTom heat load 
the SSDR-2 will calculate, as an option, the minimum possible heat load 
for the optimum :system design for each of the control methods. 

This routine option jls useful for aesign problems in which the system minimum 
heat load is not defined. The control methods are, in many instances, 
integrated with panel flow patterns and as shown in Table 3-2 are 
generally in order of increasing heat load range and control method com- ; 
plexity. The definitions . of the seven Table 3-2 headings are presented 
in Table 3-3. 

The two-dimensional panel characterization used in. SSDR is 
shown in Figure 3-2. The method of characterization assumes a triangular 
panel for panel designs where the required area is less than. one-half the 
available area. Two triangular panels are assumed for required area greater 
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TABLE 3- 2 


PANEL /SYSTEM HEAT LOAD CONTROL 
(See Table 3-3 For Definition of Headings) 



Panel Flow 
System 

H 

eat Load Control 

Configuration 

Series 

Parallel 

Bypass 

Regen. 

Valve 

Stag. 

1-D Tube 
Pattern 

2-D Tube 
Pattern 

Notes 

1. 

X 


X 



X 



2 

X 



X 


X 



3 

X 


X 

X 


X 


1 

4 


X 

X 



X 



5 


X 


X 


X 



6 


X 

X 

X 


X 


1 

7 


X 



X 

X 



8 


X 



X 


, X 

2 

9 . 


X 

X 


X 

X 



10 


X 

X 

i 

X 


X 

2 

11 


X 


1 

1 

X 

X 

X 



12 

■ 1 

X 


X i 

X 


X 

2 

13 


X 

X 

1 

X 

X 

X 


3 

14 


X 

X 

X 

X 

j 

X 

3,2 


1. The Bypas s -Regenerative control method will be considered only if 
Bypass control and Regenerative control methods are not feasible. 

2. Two dimensional tube pattern radiator panels will be considered only 
if one dimensional tube pattern panels are not feasible for this control 
method. 

3. The Bypass -Regenerative -Stagnation control method will be considered 
only if Bypass -Stagnation control and Regenerative-Stagnation control 
methods are not feasible. 
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TABLE 3-3 

HEAT LOAD CONTROL AND PANEL FLOW PATTERNS 


PARALLEL/SERIES 
FLOW PATTERNS 




SEno seeeentIne plow 


Parallel Flow Paths may be joined by large 
low velocity headers or central flow distribu- 
tion manifolds depending upon the type of fluid 
used and the method of heat load control. 




BYPASS AND 
regenerative HEAT 
LOAD CONTROL 



BYPASS CONTROt 


At low load, fluid is bypassed around the 
panel to reduce average radiator temperature, 
thereby reducing the heat rejection. The 
design limit is reached when the fluid outlet 
temperature hears its freezing point. 



REGENERATIve CONTROL 


At low load, the temperature of the fluid at 
the radiator inlet is reduced to achieve a lower 
average radiator temperature as the fluid outlet 
temperature nears its freezing point at the 
design limit. This provides wider heat load 
controT than bypass control but adds another 
system component which can become quite 
heavy. System pressure drop with full flow 
through the radiator at low load also limits this 
design. 



COMEINEO BYPASS AMD RECCMBRATIVI CONTBOL. 


A combination of regeneration and fluid bypass 
with two valves may be utilized to decrease the 
size of the regenerator a,nd improve pressure 
loss at low load, but logic to transfer control 
from the bypass valve to the regenerator flow 
control valve must be provided in the system. 
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TABLE 3 -3 (Continued) 

HEAT LOAD CONTROL AND PANEL FLOW PATTERNS 


STAGNATION HEAT 
LOAD CONTROL 


HAPtAtor 



VALVe STACHATIOM 


An automatic or manually operated valve can 
be used to shut-off flow in a bank of tubes to 
decrease radiating fin effectiveness. The 
design may allow for fluid freezing in the 
non-flowing tubes if necessary to achieve a 
sufficiently wide variation of radiator heat 
rejection. The warm tube(s) on the panel can 
be used to proyide heat to the frozen tubes for 
resumption of full flow. This valve stagnation 
method can be used in addition to system bypass 
and regenerative heat load control. 


EFFECT OF TUBE 
PAT TERN . ON HEAT 
LOAD CONTROL 



OME OIMEMSIONAL RADlATINC FIN 


Parallel tubes which flow in one direction 
oh a rectangular panel yield a low load 
stagnated flow temperature profile which 
can be , approximated by one dimensional fin 
theory to calculate low load performance using 
L as the fin half width. A limiting low load 
fin effectiveness is reached as a result of 
limitations on panel thickness and the speed with 
which full flow for high load heat rejection must 
be re-established- 


FLOWING tUBE 
AT LOW LOAD 



TWO DIMENSIONAL RADIATING FIN 


Parallel tubes in which flow turns at 90° on a 
panel will yield a deeper low . load stagnated flow 
temperature profile which can be approxirhated 
by radial fin theory. For siniilar overall panel 
geometries the radial fin will provide twice the 
heat load, control range. These tube patterns 
can be used with any of the methods described 
above. 
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than one-half the available area as shown in Figure 3-2. With this procedure, 
the tube spacing remains constant for the entire radiator panel. 


3-2.2 Description of Calciilation Methods 

A brief description of the analytical techniques and calculation 
methods are presented in this section. 


The internal convective heat transfer coefficient is calculated 
assuming fully developed flow conditions for the two flow regimes: transi- 

tion and turbulent. For the laminar flow regimes an average value of heat 
transfer coefficient is used which includes that for the entry length region 
of flow for each panel. Assuming constant temperature per unit length, 
the convection coefficient is given by: 

h = I [Nu] 

where : 


h = the internal flow convection coefficient 
k = . fluid thermal conductivity 
Nu = Nusselts Number 

D = the hydraulic diameter of the flow tube 


The commonly used equations for determining the turbiilent and 
transition flow Nusselt Numbers are used in the computer routines. For 
the laminar flow regime a curve fit has been made for the mean Nusselt 
Number which includes the effect of entry length (Ref 2 ), . The fluid pro- 
perties used are the . effective fluid properties as defined in Section 3.2.3- 

For Reynolds Nimibers less than 1900 the curve fit for the 
mean Nusselt Numbers used: 


Nu = 3.65 + 


.0668 (D/X) RePr , 

1 + .0V[(D7x)RePr ] 2/3 


For Reynolds Numbers greater than 6 U 00 
equation for turbulent flow is used: . 

Nu = 0.023 R 

e 


the Dittus-Boelter 


For Reynolds Numbers between I 9 OO and 6 UOO, a modified H. Hausen 
expression for flow in the transition region is used: 

Nu = 0.166[R^^'^^ - 125 ] Pr^'^^ 

where : . 

Nu = Nusselts Number 

^ M ’ • 

: k 
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X = the distance from flow entrance 

Re * ReyiioldB NVimber 
Um. 

~ IIP . 

m = . mass flow rate 

ft = viscosity of the fluid . 

P = wetted perimeter 
Pf = Prandtl Number 


Gp = specific heat of the fluid ' 

In many instances, panels with stagnation control techniques 
require the prime flow tube (that radiator tube with the. highest flow during 
low load stagnation) to have a larger diameter than the other tubes for optimum 
design. This occurs because pressure drop constraints at low load often 
requires a larger diameter than the optimum for the remainder of the panel. 

In SSDR the. user has the option to specify either thkt all tubes be identical 
or that all tubes except the prime tube are identical. 

The radiator- panel pressure loss is calcvilated by the Fanning 
equation for a single phase Incompressible fluid. The fluid properties 
used iu the pressvure dTpp equiatibh are the affective fl\iid properties as 
defined in Section 3.2.3. The pressure drop for each panel is calculated 

AP = f L/D p V^/2g 


relation: 


The SSDR calculates the friction factor for laminar flow by the 




For turbvilent flow the friction factor is calculated by;: 


For the transition region a pplynominal curve ;fit; is used to remove any 
discontinuities from the friction factor Ciirye: 

. f A Re + B Re + C Re + D 



The constants. A, B, C 
(1) f = f 


(2) 

(3) 


f = f 
dRe 


(M 


df 

dRe 


L 

t 


and D can he determined from the four conditions : 
at Re = Re,, (beginning of the transition region) 
at Re = Re. (end of the transition region) 

Xt 

-64 (Re^ )”^ at Re = Re, 


-.0T9 


at Re ” R®.|. 


where : 


f = friction factor 

L = panel flow length 

D = tube diameter 

Pe = effective fluid density 

V = fluid velocity 

Re = Reynolds N\unber 

Re, = laminar end of transition region 

Re. = turbulent end of transition region 


The radiator fin effectiveness is defined as the ratio of the 
heat rejected by the actual fin to the heat which wo^lld be rejected if the 
fin were at the fin base temperature. Fin effectiveness is dependent on 
the fin geometry (length, thickness), the fin material properties, the 
base temperature, and thermal environment. Values of fin effectiveness 
have been obtained (and presented ih the literature) by n\mierically 
integrating the non-linear fin heat flux differential equations. Curve fits 
of fin effectiveness for flat rectangxilar fins were obtained, in Reference 
(3) and are employed for the SSDR effectiveness calculations. 

The characterization of a heat absorbing fin is included in the 
SSDR. The differential equation describing the heat flow in a ohe-dimensional 
radiating fin was integrated for a reasonable range of boundary conditions 
to obtain effectiveness for values of sink temperature greater than the fin 
base temperature. A ci^ve fit of the information, was performed and the 
results programmed into SSDR. An equation was derived that results in 
effectiveness values with a +1^ degree of accuracy: for sink to base tempera- 
ture ratios under a prescribed boundary. For a detail discussion of heat 
absorbing fin characterization see Appendix B. , , 
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The rectangular fin effectiveness equations are a function of 
the sink-to-base temperature ratio and the rectangular fin parameter: 


T 


= (T /T.) 
s b 


)| 


2 3 

Av = ^ 


KZ . N 


surf . 


where: . 

= sink to base temperature ratio 
Xv = rectangular fin parameter 

Y = fin length 

^ € = thermal emissivity .■ / 

a .= Stefan - Boltzmann constant 

= fin base temperature, °R 
K = thermal conductivity, BTU/hr ft^R 

Z = fin thickness, ft 

N^urf ” number of fin radiating surfaces (l or 2) 
= effective sink temperatiire , °R 


Circular fin effectiveness theory is used for radiator panels 
operating at minimum heat load and employing two dimensional tube, pattern - 
stagnation heat load control methods (see Table 3-3). Values of circular 
fin effectiveness were obtained by numerically integrating the non-linear 
heat flux differential equations. Curve fits were made and included in SSDR 
effectiveness calculations. Circular fin effectiveness determination is 
presented in detail in Appendix A. 


The fin effectiveness equations: for circular fins are a function 
of the fin parameter and the radius ratio: 


Ac =, (rp - 


KZ N 


sur f ■ 


R = ^o 
r • — 

r 

1 



where : 


Ac = circular fin parameter 

R = radius ratio 
r 

r^ = radius of fin outer edge 
= radius of fin inner edge 

Micrometeoroid protection req^uirements for the radiator panel 
tubes are considered in radiator weight calculations by the use of a tube 
thickness equation. The tube thickness is a function of the probability of 
penetration, the exposed tube area, and the mission duration. The thickness 
equation will thus calculate the tube wall thickness necessary for radiator 
system compatibility with the mission definition. The user may specify a 
minimum tube thickness (which overrides the calcTilated thickness) to be . 
considered in the event that the micrometeoroid protection requirement is . 
not as severe as the required panel manufacturing constraints. The equation 
for the calculation of the required tube thickness is given below: 

WPTP 

t. = WPTl X (Time x A.) (3.1) 

U . u ■ ■ 

where: 

t . =. tube thickness - ft 

Time = mission duration - days 

■ 2 

At ' . = total tube projected area - ft 

WPTl, WPT2 = equation constants 

The radiator area required to reject the design maximxim heat 
load is calculated by iterating on the area and heat balance between the 
fluid, tube, and panel until the required heat rejection and temperature 
performance is obtained. The base temperatiire. is calculated, by equating the 
heat radiated by the panel to heat lost from the fluid through ah iterative 
process . For the conditions where the required heat rejection cannot be met 
due to area limitations, the excess heat not rejected by the system is assumed 
to be rejected by a water boiler/sublimator . The heat balance equations 
used in the iterative process is given below: 

Panel - (^ - ij A[ '^T^^ Q^] 

Fluid - = (i) C- (T. - T ) 

. . P out 

Tube Op = hA(T^ - Tg) , 


where: 


Qp = heat rejected 

T] = fin effectiveness 

( - thermal emissivity 

o - Stefan - Boltzman constant 

T_ = base temperature 

= heat absorbed, by panel 
ti> = system flow rate 


T. 

in 


out 

hA 


= fluid. spec if ic^heat 
= panel inlet temperature . 
- panel outlet temperature 

= heat transfer coefficient 

T, +.T , 

= . in out 


For control methods where a regenerator . is used. the heat rejected by. the 
regenerator is calculated as shown below. 

SeGN ~ ■ ^out^ 

3 . 2.3 Determination of Fluid and Metal Properties 

Fluid Properties 
Thermal Conductivity 

The fluid thermal conductivity is calculated by: 

K^(T) = KFO + KFl X T + KF2 x T^ + KF3' x T^. (Equation 3.2) 

where: 

K„(T) = thermal conductivity r ^77 5 ^ . 

.1 n2* X 1 "C X K 

T =. temperature - °R 

KFO, KFl, KF2, KF3 - equation constants 



where : 


where : 


The effective fluid thermal conductivity for a radiator panel is 

° ^ T. 

1 

Density 

The fluid density is calculated by: 
pf(T) = DFO + DFl X T + DF2 x + DF3 x (3.3) 


pf(T) = fluid density - lb /ft^ 

m 

T . = temperature - °R 

DFO, DFl, DF2, DF3 equation constants 

The effective fluid density for a radiator panel is 
■ T 

"o-^l J 

i 


P^(T) dT 


Specific Heat 

The' fluid specific heat is calculated by: 

C (T) = CPO + CPI X T + CP2 X + CP3 x 

P , 4 ■ 1+ 

■'+ -CP X T • 


(3.4) 


e (T) = specific heat - BTU/lh 
P ^ m ' 

T . = temperatiire - °R 

CPO, CPI, CP2, CP3» CPU - equation constants 


The effective specific heat for a radiator, panel is : 
C 

pe T 


/ 

o 1 • 


^ C (T) dT 
P 


T. 

1 • 


Dynamic Viscosity 

The fluid viscosity is calculated by: 
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(3.5) 


1At)> 

where: 


(MUO+MUl/T + ^ ^ 

'r 



M (T) = fluid viscosity - Ib^/ft-hr 


T ' = temperature - °R 

MUO, MUl, MU2 - eciuation constants 

The effective fluid viscosity for a radiator panel is; 


T . 

(T) dT ■ 

Metal Properties 

The metal conductivity is calculated hy: 

K (T) = KMO + KMi X ,T + ,KM2 X T^ + KM3 X T^ (3.6) 



where: 

K (T) =thennal conductivity - BTU/hr-ft-®R 
. m 

T . = temperature - °R 

KMO, KMl, KM2, KM3 - equation constants 
The metal density is calculated by: • 

p (T) = DM0 + DM1 X T + DM2 x T^ + DM3 x T^^ (3-7) 

where : 

P (T) - metal density - lb. /ft^ 

■ - m . ' ' m ■ • . 

T = temperature - °R 


DM0, DM1, DM2, DM3 - equation constants 
3. 2. It Weight Penalty Equations 


The weight associated vith a radiator or absorber system consists of 
the tubes, manifolds, interconnecting tubing, fluid, radiating fins, control 
valves, regenerator, panel Support structure, weight due to pimping 
power, and. excess heat penalty diie to water boiling requirements at high 
heat load. The equations, used in SSDRrr2 to Calculate radiator system weight are 
presented in this section. . 



Radiator Tube WeiRht : 

W^, = N X X L x ff X t^, (D^ + t^) (RHimiNP + 1) , (j. 8 ) 


Manifold Weight : 

W„ = NMMT X P ,, X Y X X t,(D„ + t.J (REDUNP + l) ( 3 . 9 ) 

M M MM M 

Interconnecting Tube Weight : 

X XLL X 7T X tj^ (D^ + t^) (REDUNP +l) (3.10) 

Tube Fluid Weight : 

= : N X Tg X L X n X (REDUNP + l) (3.11) 

Manifold Fluid Weight : 

Wj^ = NMANT X ~p X J X n X (REDUNP + l) (3.12) 

r 

ConnectinR Tube Fluid Weight : 

='^ X XLL X 7 T X D^^ . (REDUNP + 1) (3.13) 

Iir » e • • L 

T ■ 


w„ 


= weight of tubes per panel - lb 


m 


W, 


M 




W, 


TF 


W, 


MF 


= weight of manifolds per panel - lb 


m 


= weight of connecting lines - lb 


m 


= weight of fluid in tubes - lb 


LF 

N 

NMANT 

Pe 


m 

weight of fluid in manifolds - lb 

m 

weight of fluid in connecting lines 

number of tubes of panel 

number of manifolds per panel 

density of metal - lb /ft 

m 

effective fluid density - lb /ft 

m 


- lb 


m 


L = length of tubes - ft 

Y = width of panel - ft 

XLL . = connecting line length upstream of panel - ft 

t^ = tube thickness - ft 

t., - manifold thickness - ft 

t_ = connecting line thickness - ft 

. .. L 

= tube inside diameter - ft . 

D-, = manifold inside diameter - ft 

D = connecting line_ inside diameter - ft_ _ _ , 

... jj . .... . ^ _ . _ _ . 

REDUNP = number of redundant fluid systems 

Radiating Fin Weight : (See Figure 3-3) 

W„: - P X X X Z X Y 

fin m 

= weight of radiating fin - Ib^ 

■ 3 • 

p = density of metal - lb /ft 

X = length of fin - ft 

Z = fin thickness - ft 

Y .= width of fin - ft, [Y panel - ,N(D. + 2 t, )] 
Control Valve Weight : 

■ = NVT X WPV X (REDUNP + l) . 

where: 

W = weight of valves per panel - lb 

V , , m 

NVT = number of valves per panel 

WPV = weight per valve - Ib^ 

REDUNP — number of redundant fluid systems 
Regenerator Weight ; 

Wj^ = X WPR 


(3.14), 


(3.15) 


(3.16) 
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where : 


where : 


where : 


where : 


= regenerator weight - lb 
R . . m 

= heat transfer across regenerator - BTU/hr 
WPR = regenerator weight penalty constant - lb^/(BTU/hr ) 
Radiator Structiiral Support Weight : 

= WPSl + WPS2 X A + WPS3 x (3.17) 

Q 

= structural weight - lb 
S m 

2 

A _= panel area - ft 

WPSl, WPS2, WPS3 - structixral weight penalty constants 
Pump Power Weight Penalty : 

W„ = WPPP X Af yi 4./ ~ X ll+lt (3.l8) 

P tot e . 

Wp ' T puinp power weight - Ib^ . 

= pressure drop - psi 

W . , = total coolant weight flow - lb /hr 

WPPP = pump power weight penalty - lb^/(ft-lb^/hr V 

_ = effective coolant density - lb /ft^ 

p e m . 

Excess Heat Rejection Weight Penalty : 

= WPEH X Time X ((^ - C!^ej> (3.19) 



= excess heat rejection weight - 

WPEH = weight penalty - Ih^/ (BTU/hr) - days 

Time = mission duration - days 

= required heat rejection - BTU/hr 

Q . = actual heat rejection - BTU/hr 

rej “ ' . 
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ROUTINE OPERATION 


The calculation procedures and the overall operational flow of 
the Steady State Design Routine are presented in this portion of the report. 

The SSDR can be conveniently partitioned into ten major sections 
which are functionally independent. These major sections are defined as: 

(l) input, (2) initialization calculations, (3) panel geometry 
estimate, (U) system area estimate, (5) low heat load feasibility 
determination, (6)- maximum heat load performance and area determination, 

(7) weight calculations, (8) system low load performance calculations, 

(9) optimum design and off design parameter storage, and (lO) output. 

The SSDR program operational flow chart, showing the relationship between 
the major sections, is presented in Figure 3-^* The calculation procedures 
and operations of each major, program section shown in Figure 3-^ are sum- 
marized below. 

(lOO)* Input - Reads the data input from cards, checks for input errors, 

and writes diagnostic messages for any errors. The data, card input required 
is described in' Section 3*^«^5 and a resume of input errors (fatal and non- 
fatal) is presented in Section 3.^.5* . 

(200) Initialization .Qalcul-ations. ^ - Limits for the parameters varied 

during radiator system optimization are set, sink temperatures and average 
environments are calculated, ..and the input units are converted to the FPH 
system. The fluid properties curve coefficients are called from the pro- 
perties library or from the input data. The flow rate at maximimi heat load 
and the system inlet temperature at minimum heat load are calculated. The 
limits of radiator panel fluid temperature and flow rate at minimiim heat load 
are also set for each of the various control methods. • 

(300) Panel Geometry Estimate - The tube diameter, number of tubes, and 

the number of panel flow paths are estimated for each panel type for use in 
the system are estimate. The average of . the parameter limits is used to 
determine the diameter and the number of tubes. If either or both of these 
variables has been specified in the input, the specified input value is 
used. The number of flow paths is set equal to the number of tubes (parallel 
flow). The pressure drop is calculated using the chosen geometry and the 
total radiator area available. If the panel pressure drop design requirement 
is exceeded, the n\amber of tubes and/or the tube diameter is increased until 
the pressure drop requirement is satisfied. 

(UOO) System Area Estimate - The radiator system area , is estimated 

for use in determining the feasibility of the various low heat load control 
methods. The above panel geometry, the average thermal environment , an 
average radiator panel fin effectiveness, the fluid to tube. heat transfer, 
and the limits of panel fin thickness are used to iteratively calculate the 
radiator panel area. 

* Numbers refer to the statement numbers in the Program and on Figure 3-^. 
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Repeat for 
Additional 
Fluids and 
Metals 

~s~ 


100 * 


© 

faput [ 


'400 



. 

llhitialization Calculations I 

— ^ 

■ 


^Numbers refer to 
Statement numbers in 
The Program 


300 


IPanel Geometry Estimate | 


400 


System Area Estimat^ 


500 


I 


Low Heat Load Control Method 
Feasibility Determination 



J 

High Heat Load Performance 
Arid Arigji Determination 


700 


j^Weight Calculations j 


Variables 

Incrimented 

1, tube diameter 

2. fill thickness 

3i number of tiilics 

4. number of flow paths 

5. control method 



840 


.. met j 

r__ 


Optimum or Off Design Parameter Storage 



t » 


880 


Output of Design and Off Design Data 



FIGURE 3-1+ SSDR PROGRAM FLOW CHART 
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(500) Low Heat Load Control Method Feasibility Determination - The 

feasibility of each of the eleven low heat load control methods is determined 
unlnf.’i the estimated nystem area - geometry, the control method tempei'uture 
and flow limits, and the minimum thermal environment. The radiator system 
performance must meet the minim\im heat load requirement and pressure drop 
requirement if the control method is to be considered feasible and is to 
be used during the weight optimization procedures that follow. Unless 
specified by input, the combination control methods (i.e. bypass - regenerative) 
will not be considered if the less complex control methods are feasible 
(see Table 3-2). 

The subsequent sections of the program begin the radiator system 
optimization process. As shown in the flow chart of Figure 3-^5 the 
weight optimization is accomplished by: calculating the required radiator 

area and performance at maximum heat load for a given set of geometry para- 
meters and heat load control method (6d0); determining the system weight 
(too); and, if the system is optimum or within the. specified percentage of 
the previous optimimi, calculating the minimum heat load performance (735)- 
If the system then meets all the performance requirements, the design is 
stored within the program for later printout (8Uo). The optimization process 
above is performed parametrically by systematically varying the panel geometry 
parameters of tube diameter,, fin thickness , number of tubes (tube spacing), 
and the number of flow paths. Details of these five program sections are 
described below. 

(600) High Heat Load Performance and Area Determination - The radiator 

area required to reject the design maxim\im heat load is calculated (in this 
section) by iterating on the area and heat balance between the fluid, tube, 
and panel until the required heat rejection and temperature performance is 
obtained. If the pressure drop limitation is not met in the performance 
calculations, the parameter variables are incremented and the procedures 
are repeated again; For the conditions where the required heat rejection 
cannot be met due to area limitations , the excess heat not rejected by the 
system is assTimed to be rejected by a water boiler/sublimator . 

(too) Weight Calculations - The system weight due to the radiator 

panel fins, tubes, fluid, manifold, interconnecting tubing, control valves, 
supporting structure, pump power weight penalty, and excess heat weight 
penalty is calculated using the equations of Section 3.2.^. 

The ratio of heat rejection to weight is calculated. If this 
ratio is greater than the optimum or within a specified percentage of. 
the optimiom (off design condition) , the program proceeds with the low heat 
load performance calculations. If the ratio is not within the specified 
percentage of the optimum, the parameter variables are incremented and the 
routine returns to the high heat load section for additional calculations. 

(T35) Low Load Performance - The minimum heat load radiator perfor- 

mance is determined by iterating on radiator panel fluid temperatures or 
panel flow rate (depending on the heat load control method) until the heat 
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balance between the fluid, tube, and radiating panel is obtained. The 
calculated performance must meet the design heat rejection requirement, 
within the minimum fluid temperature limitation, and must also, meet the 
panel pressure drop requirements if the radiator system is to be considered, 
an optimum or an "off" design system. If these requirements are not, met, 
the parameter variables are incremented and the routine returns to the high, 
heat load calc\ilations . 

For the routine option which calculates the minimum heat load 
possible by the system, the program iterates on fluid temperature or flow 
rate, and heat load until the heat balances are obtained. The pressure 
drop requirement at the minimum heat load is not considered, for the minimum 
heat load determination option. 

(84o) Optimum or .Off Design Parameter Storage - After the optimum 

or off design radiator system has been configured and has met the low heat 
load performance requirements, the area requirements, geometry parameters 
( fin thickness , . tube, diameter , ,etc . ). , _and majcimum and minimum„heat load 
performance predictions are .stored in data arrays (for each control method 
and panel type) for later output. If a new optimum system is calculated 
during the systematic parameter variation, the. old optimum eind each of the 
off design systems are checked to make sure that they are within the specified 
percentage of the hew system heat . to weight ra,tio. . If they are noti the old 
designs are purged from the. off .design storage arrays. Up to fifty off 
designs can be stored' in the routine. 

(880) Output of Design: and Off Design Data The output of the optimum 

radiator system and the off design systems'-.fbr .each control ' method is accom- 
plished in this part of .the routine. 

3.1+ ' USER'S MANUAL 

The Steady State Design Routine' (SSDR).- was written in FORTRAN V 
for use on the NASA-MSC Univac II08 computer. A description of the Uni'vac 
1108 computer requirements , an explanation of the input data required, and 
a SSDR run failure analysis are presented in the sections that follow. 

3.I+.I Univac II08 Subroutine Requirements 

The following is a list of the Univac ,1108 , .FORTRAN V, system 
subroutines which are. required for use with. the SSDR: , 


1. 

NERR$ ■ . 

. 13. 

. DEPTH 

2. 

■ NEXP6$ 

. II+ . 

' nfmt$ 

3. 

NXPAF$ • . 

15. 

NIER$ 

>1. 

NypAxI ; ; 

. If. 

NTNPT.t 

» 1- 

li!'XT;.'. , 


||'l.,UA'.rx 

*6. 

Ai..00 ' 

18. 

NEXP$X 

. T. 

N0UT$ 

19. 

, CONVTX 

8. 

■ NTAB$ • 

20. 

' NININ$ 

9. 

NFTV$ 

*21. 

SORT 

10. 

NI0lN$ 

*22. 

TANK. 

11. 

NOTINi 

23. 

NST0P$ 

12. 

FPACK$ 




* These five subroutines are necessary regardless of the system 
on which the program is run. 
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3.^.2 MSC Run Submission Requirements 

For operation on the MSC Univac 1108 system (FORTRAN V), the 
SSDR program is stored on magnetic tape and the data deck, with the appropriate 
monitor controls are submitted on punched cards. The order of punched card 
data deck and control card set up is as follows : 

$ JOB 

1 

Q - ASG - A = XXXXX (Progreim Tape Number) 

J - XQT - CUR 

— TRW-A 
— IN-A 

— -TRI-A .... 

Q - XQT - DECK 1 

DATA CARDS (See Section 3.4.4) 

8 - EOF 

is the tape on which the SSDR program is stored and is always 
an input tape. (The number of the ”A” tape should be punched immediately 
following the equal sign without skipping any spaces.) 

The data deck and control card set up must be accompanied by 
a MSC Form 588 run request card. The input tape number and the problem run 
time required are designated on this run request card. 

3.^.3 Run Time and Output Estimation 

The run time for the SSDR may be estimated for the Univac 1108 
by using the following equation: 

RTIME = 0.6 X NFLUID x NMETAL x [DMA^DMINj 

r ZMAX-ZMIN , 1 
DZ / ^ STMIN 

where : 

RTIME - computer run time, minutes 
NFLUID = number of fluids considered 



NMiSTAIi = number of metnla coiiBldered 

DMAX = maximum tube diameter considered, in ' 

DMIN = minimum tube diameter considered, in 

DD = diameter increment, in 

ZMAX = maximum panel thickness considered, in ^ 

ZMIN -= minimum panel thickness considered, in . , 

DZ = panel thickness increment, in , 

STMIN = minimum tube spacing, in 

, - For the cases using the '--maximum and minimum -values -for tube 

diameter, panel thickness, and tube space stored within the SSDR, the run 
time may be estimated from: 

RTIME = 10.0 X NFLUID X NMETAL 

The expressions for run time estimation are approximations since 
the amount of time spent during, convergence for high heat load panel optimi- 
zation and on low heat load performance calculation may vary considerably from 
design to design. . . 

The printed output from the SSDR may be, estimated by 
NPAGES = T X NOON x NFLUID x NMETAL x NPT +3 

where: 

NPAGES = number of pages of printed output 

NOON = number of low heat load control methods to be 

considered 

NFLUID = number of fluids to be considered 

NMETAL ■ = number of metals to be considered 

NPT = .number of panel types to be considered 

. The number of pages of printed output may be somewhat less 
than estimated from' the above eq,uation depending on how many designs are 
within the specified percentage of the weight of the optimim design. 
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3 . ^ ^ SSDR Data Card Preparation 

An explanation of the input data required by the SSDR is 
presented in this section of the report. The flow chart of Figure 3-5, 
showing the data card order and requirements, is intended to supplement 
the input description of the following petragraphs. 

Floating point data input is designated by the format specifica- 
tion ”F". The decimal point for this data may be written in any column of 
the field and its position will override the indicated position in the 
format . specification. Interger data input, as designated by the Format 
"I”, must be right adjusted within the specified field. 


Columns 

Format 

Fortran 

Nomenclature 

Description 

Card 1 

(Title Card) 



1-75 

15A5 

TITLE 

Any 75 alphameric characters to 
be used for page heading. 

Card 2 




1-10 

FIO.O 

QMAX 

Design maximum heat load, BTU/hr. 
QMAX > 0 for radiator systems. 

QMAX 4 0 for solar absorber systems. 

11-20 

FIO.O 

QMIN 

Design minimum heat load, BTU/hr. 

If blank, program will calculate the 
minimum heat load possible. 

QMIN > 0 for radiator systems. 

QMIN 4 0 for solar absorber systems . 

21-30 

FIO.O 

TINF 

Radiator system inlet temperature 
at design maximum heat load, °F. 

31-i+O 

FIO.O 

TMIXF 

Radiator system mixed outlet temp- 
erature, °F. 

i+1-50 

FIO.O 

DPMAX 

Maximum allowable radiator system 
pressure drop, PSI . If blank,- the 
system pressure drop constraint is 
not considered . 

51-60 

FlO.O 

TIME 

Mission duration, days. Used in 
tube thickness Equation (3.I) and 
excess heat Equation (3.12). 

61-70 

FIO.O . 

RET 

Reynolds number at which flow becomes 
turbxilent (pressure drop equations). 
If blank, program sets RET = UOOO. 
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Coliimns 

Format 

Fortran 

Nomenclature 

Description 

71-80 . 

FIO.O 

REL 

Reynolds number at which flow goes 
from bho iwnlnar to (.lio trurisltion 
region (pressure drop equations) . 
If blank, program set REL = 2000. 

Card 3 

(Constants , 

Weight Penalty 

Equations) 

1-10 

FIO.O 

WPPP 

Constant, pxamping power weight 
penalty Equation (3.11)5 
lhm/(ft-lhf/hr) 

If WPPP ^ 0, the program sets 

WPPP = 0.005. 

11-20 

FIO.O 

WPR 

Constant, regenerator weight 
penalty Equation ( 3 .9) 5 Ibj^/BTU/hr ) . 
If WPR 0, the program sets 

WPR = 0.002. 

21-30 

FIO.O 

WPTl 

Constants, tube thickness Equation 

(3.1). 

31-JiO 

. FIO.O 

WPT2 

If WPTl L. 0, the program sets 
WPTl = 9.16 X 10"^, and WPT2 = 

0.351. 

Ul-50 . 

FIO.O 

TTMIN 

Minimum tube thickness to be con- 
sidered, inches. 

Overrides value . from the tube 
thickness equation if less tha:n TTMIN 
If TTMIN C 0, the progr^ assumes 
TTMIN = 0.035. 

51-60 

FIO.O 

WPEH 

Constant, excess heat rejection 


weight penalty Equation .( 3.12) , 
lbia/(BTU/hr)/Day, , 

If WPEH ^ 0, the program assumes 
WPEH =0.008. 

If WPEH = 0, the program assumes 
WPEH = oo . . . \ 

61-70 FIO.O ; WPV . .. Constat, valve weight penalty 

Equation (3-8), lb / valve. If 
WPV 4 0, the program sets WPV = 1.0. 

T1-80 FIO.O REDUNP Number of redundant ‘ fluid systems . 



Fortran 

Columns Format Nomenclature Description 

Cord h (Constants, Weight Penalty Equations) 

1-10 FIO.O ,T0L0D Weight percent from, optimum considered 

program printout off designs. 

If TOLOD ^ 0, the program sets TOLOD = 10 


11-20 

FIO.O 

WPS! (1) 

Constants , structural weight penalty 
Equation (3.8), panels of type 1. 

21-30 

FIO.O 

WPS2 (1) 

If WPS! (l) Z. 0, the program sets 
WPSl (1) =2.0, 

31-1+0 

FIO.O 

WPS3 (1). . 

WPS2 (1) = 0.2, . 

WPS3 (1) = 0. 

Ul-50 

FIO .0 

WPSl (2) 

Constants, structural weight 
penalty Equation (3.10), panels of 

51-60 

FIO.O 

WPS2’'(2) 

f ' r type 2." ' “ 

If WPSl (2) Z. 0, the program sets 

61-70 

FIO.O 

. WPS3 (2) 

WPSl (2) = 2.0, 
WPS2(2) = 0.2, 
WPS3 (2) =, 0. 

C^d 5 




1-10 

FIO.O 

DMINI 

Minimimi tube diameter to be considered 
inches . . 

If DMINI ^ 0, the program will set 
DMINI = 0.125. 

11-20 

FIO.O 

DMAXI 

Maximum tube diameter to be considered 
inches. ... 

If DMAXI ^ 0 . , the program sets 
DMAXI = 1.0. . 

21-30 

FIO.O 

DD- 

Tube diameter increment, inches. If 
DD ^ 0, the program .sets DD = 0.125. 

31-Uo 

FIO.O 

STMIN 

Minimum tube spacing to be considered, 
inches . 

If STMIN ^ 0, the program sets 
STMIN = 4.0 . 

41-50 

FIO.O 

STMAX 

Maximum tube spacing to be considered. 


inches 

if STMAX ^’0,'the program sets 
STMAX ‘ =36.0 . 


Columns 

51-60 

61-TO 

71-80 

Card 6 

I - 5 . 

6-10 

II- 15 
16-20 

21-25 

26-30 

31-35 

36-kO 


. » Fortran 

Format Nomenclature Description 

FIO.O ‘ ZfflNI Minimum fin thickness to he considered, 

inches 

If ZMINI ^ 0, the program sets 
7, MINI = O.OIO. 

I'M 0,0 '/ 1 MA.XI Maximum Vin bliirMuinss bo be eonnideretl, 

inches . 

If ZMAXI ^ 0, the program sets 
ZMAXI = 0.10. 

FIO.O DZ Fin thickness increment , inches . 

. If PZ ^ 0, the program sets 
DZ = 0.01. 


15. NP Radiator panel system index. 

NP = 1, radiator panels are in series. 

NP = 2, radiator panels are in parallel. 

15 NST (1) Number of panels of type 1. 

If NP = 2, the program sets 
; - NST (1) =2. 

’ 15 . . NST (2) ^ Number of panels of type 2. 

' If NP = 2, the program- sets 
. . NST (2) = 0. 

15 IPUMP Piamp power penalty index. 

. . IPUMP =0, penalty calculated at maximum 

, heat load 

IPUMP =1, penalty calculated a low 
heat load 

IPUMP =2, penalty calculated at average 
heat load 

15 IPRID . Prime tube diameter index.* 

IPRID > 0, prime tube diameter is calcu- 
lated 

IPRID ^ 0, all tubes .are of same diameter 

15 NFLUID Number of radiator coolant fluids 

to be. considered, Maxim\am 
limit is NFLUID = 10. 

15 NMETAL Number , of radiator metals to be 

considered. Maximum limit is 
NMETAL = 10. . 

15 IMG =0, all heat load control methods 

are to be considered. 

=1, heat load control methods 
are specified on Card 7* 
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* Colimns 

Format 

Fortran 
. Nomenclature 

Description . 

Card 7 

(Heat Load Control Method;: 

Specification)(Oinit : this cscrd for IMC = O) 



ICON (ii 

Control Index 

. = 0, do not consider . 

= 1, consider 

1-5 , 

/ 15 . 

ICON (1) 

Bypass control index 

6-10 

,15 

ICON (2) 

Regenerative control index 

11-15 

^5 

ICON (3) 

Bypass-Regenerative control index 

16-20 


ICON (1*) 
ICON (5) 

Stagnation control index , one 
dimensional tube pattern . 

21-25 

15 

Stagnation control index, two 
dimensional tube pattern 

26-30 : 

V '15; 

ICON (6) 

Bypass-Stagnation control index, 
one dimensional tube pattern 

31-35 

15 

ICON (7) 

Bypas s-St agnati on control index , 
two dimensional tube pattern 

36-1*0 

15 

ICON (8) 

Regenerative-Stagnation control 
index, one dimensional tube pattern 

1*1-1*5 

/:'■ ■:i5' 

ICON (9) 

Regenerative-Stagnation control 
index , two dimensional tube pattern 

1*6t50 

'^5 ■ 

ICON (10) 

Bypass-Regenerative-Stagnation • 
control index , one dimensional tube 
pattern 

51-55 

. ' ' ^5 . 

ICON (11) 

Bypas s -Regener at i ve -St agnati on 
control index, two dimen s i onal tube 
pattern ■ ■ ■ ' . ' 

Card 8 

(Peuiel Type: 

Data Card) 

■■•V’ ■ •; . ' 

Panel Type 1 



1-10 

FIO.O 

AMAX (l) 

Maximum , eirea available each , 
, panel , ft^ 

li-20 

FIO.O 

YP (1) < 

Panel width, ft. 

21-30 

FIO.O 

E (1) 

V Panel thermal emissivity 

31-1*0 

FIO.O ' 

FT (l) 

Panel tube emissive factor. 

(ft (1) = 1.0 for tube projected 
area as the effective relating ■ 
' area . ) If blank , program 

id. 11 set FT (1) = 1.6. 

1*1-1*5 . 

P5.0 

DS (1) V 

Specified tube diameter, inches . 
If blaiik, program ’^11 opbi^ 


.tube dimeter.' 





Fortran 


Columns 

Format 

Nomenclature 

Description 

U6-50 

F 5.0 

,ZS (1) 

Specified fin thickness, inches. 

If blank, program will optimize fin 
thickness . 


III 

NSURF (1) 

Number of panel surfaces radiating. 

55-58 

lU 

NTS (1) 

Specified number of tubes. If 
blank, program will optimize the 
number of tubes . 

59-60 

12 

NFPS ( 1 ) 

Specified number of parallel flow 
paths on the panel. 

= 0, program will optimize the 
number of flow paths 
= -1., program will consider 

all tubes connected in parallel 
= 1, program will only consider 
serpentine panel flow 
= 1 , - 1 , 0 , NITS (1) should divide 
into NTS (1) evenly. 

61-65 

F 5.0 

DLL ( 1 ) . 

Diameter of pahel fluid connecting 
lines , inches . 

66-70 

F 5.0 

. DMAN ( 1 ) 

Manifold diameter of cylindrical 
manifold, inches . 

71-75 

F 5.0 

THAN ( 1 ) 

Manifold material thickness, inches. 

76-80 

. F 5.0 

ZMINS (1) 

Fin thiclmess not attributable to 
radiator panel. (No weight penalty 
for radiator system), inches. 

Repeat ( 

Dard 8 for panel type 2 . (NST 

(2), Card 6, must be positive. ) 

Card 9 
Panel 1 
1-lt 

(Panel Data Card) 

Blank 

5-10 

A6 

PANEL (1) ' 

Panel name, any six alphameric 
characters 

11-20 

FIO.O 

QAMAX (1) 

Heat absorbed by panel from 
environment at maximum heat load, 
BTU/hr-ft .2 

21-30 

FIO . 0 

QAMIN ( 1 ) 

, Heat absorbed by panel from environ- 
ment at minimum heat load, 

BTU/hr-ft 2 . 


Columns 

Fortran 

Format Nomenclature 

Description 

3l-»»0 

FIO.O 

FI’ (1) 

Panel emissive fiicl.pf. If PP (l) 
is blank., program seta FP (l) = 1.0. 

iil-50 

FIO.O 

xll ( 1) 

Fluid connecting line length upstream 
of the panel, feet. The first pajiel 
line length should include the return 
length (down stream) from the last 
panel . 

51-55 

15 

NMAN (1) 

Number of manifolds on this panel. 

56-60 

15 

NV (1) 

Number of vadves on this panel. 

Repeat 
NST (2) 

Card 9 for 

each panel. Total 

■ntjmber of cards needed = NST (l)' + 

Card 10 

(Fluid Data Card) 


i-lo 

ilO 

iFLUlb (1) 

Index of fluid nvimber i. (See Table 
U for list of index nvanbers.) 

il-20 

FIO.O 

TMINF (1) 

Fluid 1, minimum* allowable 
temperature, °F. 

21-30 

110 

IFLUID (2) 

Index of fluid 2 . 

31-UO 

FIO.O 

' TMINF (2) 

Fluid 2 minimum* allowable 
temperature, ®F. 

Ul-50 

110 

IFLUID (3) 

Index of fluid 3. 

51-60 

F10;0 

TMINF (3) 

Fluid 3 minimvim* allowable 
temperatvLre , °F. 

61-70 

110 

IFLUID (It) 

Index of fluid k. 

71-80 

FIO.O 

TMINF (U) 

Fluid U minimtun* . allowable 
temperatxire , °F. 

Repeat 

Card iO as 

many times as needed to supply IFLUID and TMINF for each 


fluid. 


* For solar absorbers 


this is the maximum allowable fluid temperature, ®F. 


FLUID 

INDEX 

TABLE 3-4 

INDEX NUMBERS FOR THE SSDR 
FLUIDS AND METALS 

PROGRAM 

FLUID NAME • NOMENCLATURE 

POUR POINT* 
TEMP. ®F 

1 

RS89-A 

RS89-A 

-65. 

2 

MSC-198 

MSD198 

-l4o. 

3 

60/40 Propylene Glycol 

PRO-GL 

-76. 

1* 

Freon 2l 

F-21 

■ ' -211 . . 

5 

FC-43 

FC-43 

-58. 

6 

PC-75 

FC-75 

-135. 

7 

Freon E-3 . 

E-3 

-l60. 

8 

Water 

, . WATER 

32. 

9 

80/20 Methanol-Water 

ME-H20 

-65. 

METAL 

INDEX 

Metal name . , 

PROGRAM 

NOMENCLATURE 


1 Aluminum, 75S-T6 75S-T6 

2 Magnesium, AN-M-29 AK-M29 

•3 Berrylium BE 

4 Stainless Steel, SS301 ' SS-301 

5 Aluminum, . 6o6l T6 6 o61T6 
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Columns 

Format 

Fortran ,V, 

Nomenclature Description 

Cards 11 A 

, B, C, and 

D (Fluid Properties Data) 

Cards llA, IIB, IIC, and IID are 
are not, included in the program. 
IFLUID = 0 on Card 10. 

needed only for those fluids whose properties 
Cards llA, B, C, and D are needed for each 

Card llA 




l-6\ 

a6 . 

FLUID (l) . 

Fluid name , any siy alphameric 
characters ., 

Card llB 




1-10 

FIO.O 

KFO 

Curve coefficients for the fluid 
thermal conductivity Equation (3.2). 

11-20 

FIO.O 

KFl 


21-30 

FIO.O 

KF2 


31-ltO 

FIO.O 

■ KF3 


41-50 

FIO.O 

dfo ■ 

Curve coefficients, for the fluid 
density Equation (3.3). 

51-60 

FIO.O 

DFl ■ . . 


6i-T0 

FIO.O 

DF2 ' 


71-80 

FIO.O 

DF3 ■ 


Card lie 




1-10 

FIO.O 

cpo : , 

Curve coefficients for. the fluid 
specific heat Equation (3.4). 

11-20 

FIO.O 

CPl' 


21-30 

FIO.O. 

CP2 


31-40 

FIO.O 

CP3 . • 


41-50 

FIO.O 

CP4 


51-60 

FIO.O 

MUO 

Cijrve coefficients for the fluid 
dynamic viscosity Equation (3.5)* 

6i-to 

Fio.o 

MUl 


71-80 

FIO.O 

• MU2 



ho . 



Columns 

Format 

Fortran 

Nomenclature 


Description 

Card IID 






1-10 

FIO.O 

MU3 ' 

Curve coefficients for dynamic 
viscosity Equation (3.5)- 
(continued) 

11-20 

FIO.O 

mk 




21-30 

FIO.O 

MU5 




Card 12 

(Metal Data Card) 




1-5 

15 

IMETAL (1) 

. Index of metal number 
Table 3 tU for list of 
numbers . ) 

1, (See 
metal index 

5-10 

15 

IMETAL (2) 

Index 

of metal number 

2. 

11-15 

15 

IMETAL (3) 

Index 

of metal number 

3. 

16-20 

15 

IMETAL {k) 

Index 

of metal nmber 


21-25 

15 

IMETAL (5) 

Index 

of metal number 

5. 

26-30 

15 

IMETAL (6) 

Index 

of metal number 

6. 

31-35 

15 

IMETAL (T) 

Index 

of metal number 

T- 

36-40 

15 

IMETAL (8) 

Index 

of metal mamber 

8. 

hl-k^ 

,15 

IMETAL (9) 

Index 

of metal number 

9. • ■ 

U6-50 

15 

IMETAL (10) 

Index 

of metal number 

IQ. 


Cards 13A and B (Metal Properties Data) 

Cards 13A and 13B are needed only for those metals whose properties are 
not included in the program. Cards 13A and 13B are needed for each 
IMETAL = ) on Card 12. 



Columns 

Format 

Fortran 

Nomenclature 

Description 

Card 13A 




1-6 

a6 

METAL (I) 

Metal name, any six alphameric 

Card 13B 
.1-10 

FIO.O 

■ KMO 

characters - 

,G\:u:ve coefficients for the meteil 

11-20 

FIO.O 

KMl 

thermal conductivity Equation (3-6) 

21-30, 

FIO.O 

KM2 


31- UO 

FIO.O 

KM3 


41-50^ 

FIO.O ^ 

' DM0 - - ■ - - 

" Curve^coefflcients for- the metal* 

51-6o 

FIO.O 

DM1 

density Equation ( 3 -T ) - 

6i-70 

FIO.O 

DM2 


T 1-80 

FIO.O 

DM3 






3.*+. 5 


Run Failtire Analysis 


A resiame of the fatal and nonfatal errors detected during the 
data card input will he printed by the SSDR. Fatal errors will cause the 
SSDR to end further . .data processing while nonfatal errors will cause 
certain input data to be ignored. The diagnostic message printed by the 
SSDR for the fatal and/or nonfatal errors will contain the, name of the 
input variable and the card n\amber on which the variable is located. 

The fatal errors in the input data that will cause program 
termination are: 

Card 2: QMAX =0 

Card 6: NST(l) >.10 . 

NST(2) > 10 . 

NST(l). + NST(2) >10 
NFLUID >10 
NMETAL >10 

Card 8: AMAX (IT) < 0 . . 

. YP (IT) <0 
E (IT) <0 

Nonfatal errors possible in the input of data are: 

Card 10: IFLUID (l) > 9 

Card 12: IMETAL (l) > 5 

The data for each nonfatal error will be ignored by the SSDR. 

In addition to the input data errors, the routine may fail to 
design a radiator system due to design constraints imposed by the program 
user. These constraints consist. of pressure drop limitation and radiator 
area - excess heat weight penalty availability. 

If the radiator system cannot be designed at high heat load for 
the pressure drop limitation imposed (DPMAX, Card 2), the program will print 
a diagnostic message indicating the minimum pressiire drop that can be 
achieved and then terminate calculations . If the radiator system cannot 
be designed to meet the pressure drop limitation at low heat load for a 
particular control method, the program will print a diagnostic message 
indicating the minimum pressure drop that can be achieved and omit further 
calculations for the control method. 

Under certain conditions the available radiator area for the 
required high load heat rejection may not be sufficient (AMAX(IT), Card 8) 
and, therefore, a radiator system cannot be designed without the use of 
a water sublimator or boiler to reject the excess heat. If this situation 
occurs , WPEH on Card 3, of the input must be reasonable to assure meaningfiil 
system design. The program prints out of the possibility of inadequate 
radiator area exists . 
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. 4-. C) TRANSIENT PERFORMANCE ROUTINE 

This section contains the purpose, analytical ■ methods employed, 
and user's instructions for the Transient Performance Routine (TPR). 

1),1 PURPOSE 

The TPR was developed to provide rapid, easily obtainable 
evaluation of the thermal response of radiators and solar absorbers designed 
by the SSDR. The TPR results can be interpreted to provide new inputs to the 
SSDR to iterate on the initial steady state radiator design. The TPR utilizes a 
time varying environmental heat. flux simulation to determine orbital radiator 
and solar absorber performance whereas the SSDR utilizes a pair of effective 
environmental heat flvixes for minimum and maximum design points . The results 
from TPR can be used to judge the adequacy of the effective environments ini- 
tially used in SSDR; By using TPR results to adjust these environments for sub- 
sequent design iterations with. SSDR, ah optimum radiator or absorber which satis- 
fies 'transient response 'requirements “Can be evolved".” ' “ / ' 

h.2 ROUTINE DESCRIPTION AND ANALYTICAL METHODS 

The TPR utilizes a simplified LWM25 (Reference 1+) data input 
with a predetermined flow system thermal model suited to analysis of any 
number of parallel or series floirf* radiator panels . The, detailed nodal data 
needed for transient finite difference analyses of either one or two-dimen- 
sional radiator panels is generated by coupling the TPR with the Automatic 
Nodal Subdivision Subroutines (SUBANS.and SUBRAD) described in Section )|'.3.1. 

The user may easily generate nodal models of either rectangular or circular 
, radiator panels with varying degrees of fineness to establish the number of 
nodes required for accurate transient simulation of radiator perf ormanc.e . 

Any desired combination of rectangiHar and circular panels can be accommodated 
by TPR (See Figure 1-1 and 1-2) with the condition that all, parallel flow 
panels be located upstream of any series. flow panel. 

The predetermined flow system (Figure 1-1 ) provides for combina- 
tions of valves and use of a regenerator to simulate any of. the low load 
control methods described in Table 3-3. The TPR includes automatic output 
plotting, restart capability, provisions for use of cyclic incident heat 
curves and, use of data tapes and editing.. The program computational speed 
is significantly better than that of. the basic general, purpose LVYM25 
routine. A detailed description of the routine is given below. 

1.2.1 Groundrules and Assumntions 

The following. are ,the groundrules which were established by 
LTV as guidelines . for the development of the TPR. Several assumptions are 
inherently contained in the routine because’ of these grdundr^es. 

(a) Flow paths are written in to correspond. to the skeleton 
flow system (see Figure 1-1) for flow rate and inlet 
temperature. 

* All options apply to. solar absorbers, as well as radiators . 
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FIGURE 4-2 FLOW SYSTEM OPTIONS 



(b) The pressure-flow balance on each radiator panel is per- 
formed assTiming parallel flow for the panel tubes, (i.e., 
manifold losses between tubes are neglected). A minimum 
of • 2 tubes must exist on each radiator panel. Also, at 
least one tube lump and one fluid lump must exist for each 
tube . 

(c) The fluid lump number must be identical to the enclosing 
tube lump number . 

(d) There are 5 valves in the system. The valve locations in 
the flow paths are written in the routine, but the valve 
parameters are input. It is up to the user to manipulate 
the valves (by shutting off, etc.) to obtain the desired 
flow configurations. 

(e) Proportioning valve - Programming was performed to give 

three options for the proportioning valve: (l) normal 

flow proportioning, (2) no proportioning, and (3) locked 
one way for single panel operation. 


(f) The regenerator is analyzed by a steady state effectiveness 
method. The regenerator mass may be approximated by adding 
l/k of the regenerator thermal .capacitance into each inlet 
• hnbe 1 mrip and the tube 1-UTnps immediately dowhr/t;reaTn of eardi 


(g) A baseline thermal model for the flow system plumbing and 
regenerator is provided so that the user may concentrate 
on the analysis of the radiator panel design. If details 
of flow system pliambing and/or regenerator are known, or 
if the baseline model is not practical for. fluids and flow 
rates under consideration, the user can replace or alter 
the dimensions or number of nodes in the baseline model. 

U.2.2 Temperature and Flow Analysis 

Temperature Analysis - Finite difference approximations of the 
differential equations governing the temperatures for transient conditions 
are solved using the. forward differencing (explicit) method.. These equations 
and the method are described in detail in Reference 4. The assumptions that 
have been made are: 


(1) The fluid and metal thermodynamic properties are constant 
within any element for a given time increment but may vary 
between elements. 

(2) Cylindrical surfaces are approximated, by small rectangular 
segments . 

(3) Radiant interchange is not considered. 


A problem with flowing fluid is considered to have three types 
of. finite - difference lumps; i.e., fluid lumps, structure lumps, and tube 
lumps which interface with the structure and fluid lumps. The governing 
equations for each type of lump are: 


(l) Fluid Lump (The option to use an average of fluid lump 
inlet and outlet temperature with the hA .term which is 
discussed in Reference. U- is not available . in this routine) 


1 (Tfi -t/) 


+ hjAj. - Tj.") 1 


(Equation U-l) 


where: 


w. 


i 

f 

i+1 


f 

At 


w 


fu 


= weight of the fluid in Imp f 

= A^(p)(L) where p is the density, L is the lump 
length, and A^ is the Imp cross sectional area 

= specific heat of the fluid in Imp f 

= temperature of the fluid at time T 

= temperature of the fluid at time T. + 

- time increment for next calculation 

= fluid flow rate 

= temperature of fluid lump directly upstream of 
Imp f in the flow path, at time T . 

= fluid convective heat transfer coefficient 

= area for convection heat transfer 

= temperature of enclosing tube Imp at time t 


(2) Tube Imp equation 






* ? “t= <’'s" - 

+■ . (U-2) 

U8: ■ . 


T 


i+1 


tube liamp weight 
specific heat of the tube 
temperature of the tube lump at time 
temperature of the tube lump at time 'f + A 


U 


tj 


the conductance between tube lump t and adjacent 
tube lumps ' J . 


= A 


ec 


Y Y 


u 


ts 


A is the effective conduction area between tube 
cc 

lumps t and J 

Y^ is that portion of the conduction path length 
between lump t and j which lies in lump t 
Y is that portion of the conduction path length 
between lump t and J which lies in lump J 
K^. is the thermal conductivity of tube lump t 

K. is the thermal conductivity of tube lump J 

J . ’ ■ ' - 

the conductance between tube lump t» and adjacent fin 
or. structure lumps, s 


= A 


ec 


Y. Y 
- 1 + ^ 

L t s -• 



Y^ is that portion of the conduction path length 
between lump t and s which lies in lump t 
. Y is that portion of the conduction path length 
between lump t and s which lies in lump s 

= temperature of adjacent fin or structure l\imp s 
at time "t 

= incident heat on tube l\amp t at time ^ 

= tube lump external area 

= absorptivity of lump t for incident heat flux 
= emissivity of tube lump, t 


h9 


- fluid convective heat tranafer coefficient 


where: 


= area for fluid convective hea,t transfer 
= fluid lump temperature at time ^ 

( 3 ) Structure lump equation 

* :4Z[ 

wc*-, i 


"s*' 


s s 


+ a A Q ^ - f a A (T 
s s s s s s 


V ] 


(4-3) 


w 


T 


i+1 


s 

ATT 


u 


sJ 


= weight- of lump _ , _ ■ 

= specific heat of lump s 

= temperature of lump s at time *TJ 

= temperature of lump s at time'^T+ 

= time increment for next step in calculation as 
determined by convergence criteria 

= the conductance between structure liimp s and adjacent 
. structure liimps , J 


U , = A 
n J f»r; 


y. y, 

JL+.JL 

, K K, , 

L s J -1 


where ; 


T. 


Q. 


Y is that portion of the conduction path length 
between lump s and J which lies in limp s 

Yi is that portion of the conduction path length 
^ between lump s and J which, lies in Itunp J 

A is the effective conduction area between lumps . 
s and j 

is the thermal conductivity of lump's 

K l.H tiiB Lhej'mal oojiddct.l v:l l.y or lump J 

• ■ J ■ 

= temperature of adjacent tube lump J at time ^ 

= incident heat flux on lump s at time 'V 
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A = area of lump s (for radiation) 

s . . 

- absorptivity of lump s for incident heat flux 

( ^ = emissivity of lump s 

= the conductance between structure lump s and 
adjacent tube lump t 

r 1 1 



The heat transfer coefficient between the fluid and tube is 
given in Reference. U-1 as: 

Laminar Entry Length: 


h, = 1.86 

Laminar flow fully developed: 





(U-U) 


(U-5) 


Turbulent Flow: 



where: 


(»4-6) 


h^ = convective heat transfer coefficient 

= fluid conductivity 

L = length from tube entrance 

D = tube hydraulic diameter 

- fluid viscosity evaluated at fluid bulk temperature 

M = fluid viscosity evaluated at tube wall temperature 

w , 

Re = Reynolds Number 

Pr = Prandtl Number 
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= entry length heat, transfer coefficient factor (data 
input ) 

= developed flow heat transfer coefficient factor 

A value of F of Q.5T5 and F^ of 1.0 was found to fit the Graetz 
solution for constant wail, temperature, as shown in Reference U- 1 . These 
methods of calculating heat transfer coefficients were retained in Transient 
Performance Routine. . 


Flow and. Pressure Drop Analysis - The pressure drop calculations 
for the Transient Performance Routine are performed as discussed in Reference 
k- That is, the pressure drop for each fluid lump is calculated by: 


where: 


AP = 


.2 

w 


2 PA 2 
c 


,f(WP) L F 3 ( P-b/P w) 


- ,1k 


w = tube fluid flow rate 

p = fluid density 


K 


(U-T) 


A^ = fluid lump cross sectional area 


f = friction factor 

L . = 16 /Re for Reynolds niimbers less than 2000 and is read 

from input for Reynolds niambers greater than 2000 

WP = wetted perimeter 

L = fluid lump length 

= fluid viscosity evaluated at fluid bulk temperature 

= fluid viscosity evaluated at tube wall temperature 

. F^ = pressure factor for non-circular ducts; is 1.0 for 

circular ducts 

K ~ number of fluid dynamic head losses 

The pressure drop for each tube is determined by summing the 
pressure drops of the lumps in the tube. 


follows .. 


The system flow rate-pressure drop balance is determined as 
For each iteration: 
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(l) The new, flow rate entering the system is read from the 
flow rate versus time curve. 


( 2 ) : New valve positions are determined (as described in h.2.3)> 

(3) New flow rates for tubes 2 through 8 are determined as 
follows : 


Wg = I = Flow fractions determined by valve ^1 

w^ = Wg ( Xg ,= Flow fractions determined by valve #2 

^=,V2 

w^ = (1 - Xg) Wg 

■w^ •= (1 - x^) w^; 


The flow rates for w- through 1 

number of . tubes on each, radiatorf are multiplied by the 
quantity Wg /^g , • . 

new old 


(5) A check is made to see if either valve or 5 changed 
position: 


(a) If there is a change from stagnation condition to 
parallel flow condition the flow through each of 
the parallel flow tubes is Wg/n or.w^^/n depend- 
ing on whether valve 1+ or valve 5 changed positions . 

(b) If the change was from parallel flow condition to 

stagnated flow condition w. _ = w. and w._. = w... 

13 . 9 . 13+n 10 

(6) .Pressure drop. per. Irrnip is calculated by equation h-7 
and the pressure drop for all the l\imps in a tube are 
summed to obtain the pressure drop per tube. 

(7) K is calculated for each tube by 

' ^ V ' 

tube 

K — . 

W. , . 

tube 

(8) For each radiator panel in the. parallel flow condition 
(not stagnated) the flow is determined as follows: 
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(a) Calculate; 1 

1=13, 

13+n-l 

^^side 2 ■ ’'lO'^ ^ 1/K. 

i=13+n,. ^ 

13+2n-l 

(b) Calculate the new flow in each parallel path by 


w, = AP -/K, 
i panel i 


(9) Add the pressiire drop around each of the two parallel flow 
paths , i . e . , 

&Pl-= APj+ a a . - - 

AP, = AP^„. APj. 2- Ap^^* 

where: 


Ap^ = pressure drop around left, path 
,Ap^ = pressure drop around right path 
AP = pressure drop in tube 9 

Ap^q = pressure drop in tube 10 

AP^^ = pressure drop in tube 11 

Ap ^2 ~ pressure drop in tube 12 

Ap^^ = pressure drop in left side of proportioning valve 

APj^,^ = pressure drop in right side of proportioning valve 

(lO) If there is flow in both sides of the proportioning, valve, 
solve for flow around each of the two parallel paths by 
(see Proportioning Valve under Section i+.2.3)* 

. -b +\b - ac 




2 a 


where: 
a = E 


[x7 - 
•> ' VSlT 


2 E W2 
7-2— 


" > ^2 - 


ECw^r 



E = proportionality factor 


= valve position from left 
= valve position from right 
^^side 1 ^^11 


^ ^lo'^^^sid e 2^^ ^12 


w. 


10 


^10 ^2 ~ "^9 


.^11 ^9 


\2 ^10 

( 11 ) Flow rates w^^ thrpugh w^j+n-l multiplied by 


w. 


new 




old 


Flow rates ^13+2n-l multiplied by 


w. 


10 


new 


w 


10 


old 


(If flow rate is through one side only, (lO) and (11 ) are not 
performed)’ 

( 12 ) Parallel flow paths are checked for pressure drop balance. 

If all balance within a specified tolera,nce, DPTOL in data 
input, calculations are complete. 


( 13 ) If pressure drops are not balanced, the flow rate for each 
tube is averaged with previous flow rate as follows : 

w-=aw. +(l-a)w. 

1 1 1 -ij 

new old 

where: 

a is an input quantity, usually 0.5,. 


55 



(lU) Steps (6) through (13) are then repeated until a pressure 
balance is obtained. 

It. 2. 3 Component Description 

The basic model of the TPR shown in Figure H-l consists of a 
system bypass valve, a regenerator, a regenerator bypass valve, a pro- 
portioning valve, two stagnation valves, two radiator panels and connecting 
tubing. This section will describe the assumed characteristics of the 
components which make up the system. 


System Stagnation and Regenerator Bypass Valves 

For each of the two bypass and two stagnation valves (valves 
1, 2, 4, and 5) "the user has the option tOf-'specify either rate limited valves 
or polynomial valves . . 


The rate limited valves are characterized by limiting the rate 
of opening or closing by the amount shown in Figure 4-3. The set point, 
deadbadd, dX/d AT arid X shown on the figure are input values. The rate 
of opening for the valve is determined each iteration by, the difference 
between the sensor lump temperature and the set point temperature. The amount 
of valve movement for each Iteration is then the rate of opening times the 
time incremerit. 

The polynomial valve flow distribution at a btanch depends on 
the temperature of a specified lump (sensor lump) in the problem. The 
fraction bypassed is determined by a fourth order polynomial of valve posi- 
tion versus temperature. The coefficients of the polynomial are input. 

Proportioning Valve (Valve 3) - 


« • The proportioning valve (Valve 3) is designed to respond in 

th*e direction of causing the flow in the two parallel paths to exhaust at 
the same temperature. This arrangement is utilized to provide maximum 
heat rejection when two sides of a radiator system operate in a significantly 
different incident heat environment. The equation describing the . operation 
of the valve is : 

° Venous •4 t[ (’'i - Jprevlous’* ‘ ’’lt'] 

where : . ■ 


RT’ 


■LT 


Valve gain is allowed as an input 
At = time increment 
f c = valve time constant 
X^ = initial valve position 
= present valve position 
= tdmperatj^e of sensors in 
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constant in the data. 


right and left hand tubes 


X 


through Norir-Bypass Tubes 



FIGURE BYPASS VALVE OPERATION 



After the position X is determined, it is used to define the pressure drops 
in. each side of the valve through the relation: 

2 ■ 


A P = V. 

• « R.J, 


AP. 


LT 


[?] 
■■[?] 


PPARA 


E = proportionality factor ( = ^ Gp^CT 

Card D of data preparation) 


; see Parameter 


w , w_ = right and left flow rates 
nl Li 

= valve position from left 

^ valve position from right™ . ^ ^ ^ 

The valve pressure drops are. considered together with the 
pressure drops in the remainder of the right and left hand flow paths to 
determine flow rates which give a pressure balance for both sides of the 
system. Considering the pressure drop of the radiator to be a linear 
function of flow rate, such that A P = Kw , the pressure balance in the 
radiator and valve can be written as ; 

■1 2 


LXg J L^i J 


where: 

= AP of radiator right branch/right side flow rate 

= L P of radiator left branch/left side flow rate 

WpijV Wj^^ = right and left flow rates 

The pressure drop equation may be solved for the left side flow rate by 
substituting Wp^ = w^Q^ - w^^. ^ 


^ ^ ^ St .. St " StSot 


E(w„^„)‘ 

tot 

+ — 


Denoting the coefficient of w.j._ by (a), of w _ by (b), and the constant 
term by (c),,the w^^^ may be put into the standard quadratic form. 


w. 


-b +■ 


LT 


b - U a c (step 10 of the preceeding flow and 
2 a pressiure balance description) 


58 



When and are different by less, than a specified tolerance (VLVTOL, 
Card D 5 Columns 31-^»0, the value of (a) will be small and the approxima- 
tion w = c/b is used, 
hi 

Regenerator 


In. the interest of low computer run times the regenerator 
effectiveness is determined using steady state equations. If tube it is 
assumed to be side 1 and tube 2 is. side 2, the effectiveness is given by 


1 - e 


UA 
(me ) 


- t -^1 

1 L me p , 


■ (me)- . . UA 


1 _ 

(mc)^ . 


or when mc^ = mc^j 

UA 

= ■(m c .)j 


1 + 


UA 


(me), 


where: UA is the overall heat transfer coefficient 

(mc)^.is the mass flow times specific heat for side 1 
(mc)^ is the mass flow. times specific heat for side 2 

The outlet temperature on side 1 of the heat exchanger is calculated based 
on the effectiveness as follows: 

T . = T. - e (T. - T, ) 

out^ in^ in^ .in^ 

where : 

*^out ~ outlet temperature of side 1 
= inlet temperature on side 1 


= inlet temperature on side 2 

The outlet temperature on side 2 is determined by 

(l) The enthalpy on side 2 is determined by: 

h , = h. - w- /w- (h ' - h. ) 

out^ in^ 1 2 out^ in^ 
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where: 

h , = the fluid enthalpy at T . 

outg out^ 

h . = the fluid enthalpy at T . 

out^, out^ 

h. = the fluid enthalpy at T~! 

ini in^ 

w^ = fluid flow rate on side 1 

W2 " fluid flow rate on side 2 

(2) is then determined by reverse interpolation of the 

curvis of -enthalpy versus temperature for the' fluid . 

It should be stressed that this method gives the. steady state 
outlet temperatures for, the given inlet temperatures , overall heat transfer 
coefficient, and me products for sides 1 and 2. However, an approximation 
of the regenerator transient performance may be obtained by putting 1/U of 
the regenerator thermal capacitance into each inlet tube lump and the tube 
lump immediately downstream of each outlet lump. 

Radiator Panels 

The parallel flow radiator panels are each assumed to contain 
a bank of parallel tubes with one tube designated a high flow tube as shown 
in Figures i*-l and U-2. Between the high flow tube and the remaining tubes 
is a stagnation valve, either rate limited or polynomial, so that the flow 
is proportioned between the single high flow tube and the remaining tubes 
by the position of the valve. It is required that there exists at least 
2 tubes per parallel flow panel; i.e., one tube on each side of the stagna- 
tion valve. 


The TPR is setup to handle either radiators in parallel or in 
series since these are the t^es of systems which may be generated by the 
SSdS. Any combination ’of’TeWaHgular‘'T-D~,“reHaBg\ilaF'2'^D'r''Md" c“ir^^ pahels 
in a parallel or .series flow system can be analyzed by the TPR. Both parallel 
smd serpentine flow radiators can be handled with the condition that all paral- 
lel flow panels must be located upstream of any series panels. To analyze a 
series flow system containing any number of parallel and series flow radiators 
the proportioning valve (Valve No. 3) should be locked by setting the valve of 
NOP (Parameter Card c of Users Msuiual) to 1 which will force all the flow 
to side No. 2 (Figure k-2) . The parallel flow panels should be placed up- 
stream of the series flow panels and two nodes without external radiating 
area used for side 2. 
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The radiator panel nodal breakdown may be input by one of 
two ways . . 

(1) niR user may input the nodal information direoLly 

(2) The user may input overall radiator dimensions and call 

on the Automatic Nodal Subdivision Subroutines (SUBANS , SUB2D, 
and SUBRAD described in Section U.3.I) to perform the 
nodal breakdown and generate input data automatically. 

The user may also call on SUBANS, SUB2D or SUBRAD to perform the 
nodal breakdown on some panels and manually input the nodal breakdown 
for other panels. Only one run is required on the TPR to perform the 
nodal subdivision and the transient performance predictions. 

U.3 ROUTINE OPTIONS 

lt.3.1 Automatic Nodal Subdivision Subroutines (SUBANS , SUB2D, and 

SUBRAD) 

. 3 . 1 . 1 Purpose 

The creation of data decks for the LTV. Thermal Analyzer computer 
routines usually involves the tedious task of dividing a thermal model into ■ 
its component nodes , numbering these nodes, specifying their dimensions 
and neighboring node numbers , loading the resulting values on keypunching 
sheets, and assembling the decks. The models being analyzed are usually so 
varied from node to node that no systematic method can be developed by which 
they may be broken into their components; there are,, however, exceptions 
to this non-uniformity, viz., radiator panel models. Radiators usually 
involve large areas requiring many nodes to characterize them; but at the 
same time, their makeup is fairly simple, resulting in many nodes whose 
physical dimensions are the same, the only real difference being the unique 
number assigned to each node to distinguish it from the others. 

U. 3.1.2 Performance 

Three subroutines (SUBANS, SUB2D, and SUBRAD) have therefore been 
written which, when supplied with overall radiator panel dimensions and when 
given an indication of .the coarseness or fineness of the desired nodal break- 
down, will provide all information necessary to input the modelrto the computer 
routine. SUBANS is capable of i subdividing a 2-D rectangular panel, either parallel 
or serpentine floW)SUB2D hEindles 2-D rectangular panels and SUBRAD is capable of 
performing the nodal breakdown of a circular panel with parallel flow. More speci- 
fically the subroutines will divide the radiator panels into nodes, assign lump 
and type numbers to. the nodes, identify conduction and/or flow neighbors (account- 
ing for either parallel or serpentine flow), • select the proper incident heat 
curves , and provide property curve numbers — all on tape in the form of card 
images compatible with the EDIT subroutine. Thus, a user may provide nodal sub- 
division input cards for one or more radiator panels, .followed. .by a data deck 
for the remainder of the model; and in one run the radiator panels will be 
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automatically subdivided, the resulting tape will be combined with the 
input data deck, and the problem will be run to completion. 

U. 3 . 1.3 Input 

As their major input, SUBANS, SUB2D, and SUBRAD require the overall 
panel dimensions. For a rectangular panel the length in the x-dlrection (see 
Figure k-h ) and the height in the y-direction are required input and for a 
circular panel the radius of the panel is a major input item. The panel 
thickness is input for both type panels. To keep the grid size flexible 
on a rectangular panel, the user specifies the number of nodes in the x- 
direction; and in the y-direction, he gives the typical tube lump width, 
the niimber of tubes , and the number of structure lumps between tubes. For 
a circular panel the user specifies the number of tubes and has an option 
of either specifying the number of lumps in each tube or a set of limits 
for the length of all tube lumps. Starting with n lumps for the first 
(inside) tube of a circular_pan_el,_qach^ successive tube. is required to have 
some multiple of n lumps; i.e.. In, 2n, 3n, etc. Typical tube lump width 
and nmber of structure lumps between tubes are also input items for a 
circular panel. With reference to this latter number, SUBANS, SUB2D,.:and 
SUBRAD are restrained to put half as many structure lumps between the outside 
tubes and the pauel edges, thus ^the number, is- required to. be even. 

Those quantities and numbers which dp not enter into nodal sub- 
division calculations but which are required for input into the thermal 
analyzer routines are also SUBARS and SUBRAD input items; these include, 
properties and property curve numbers, initial temperature, heat transfer 
and external radiation areas, cross-sectional area, wetted perimeter, and 
friction factor coefficient. These items are simply stored until the data 
tape is being written, at which time they are output exactly as they appeared 
on the input cards; no checks are made for correctness of either the formats 
or the values. 


4.3.1*^ One Dimensional Rectangular Model Subdivision 

The performance of SUBANS is best illustrated by showin 5 
what it will do with a sample data set. Assuming a panel length, x, of 
120 inches (Figure U-5) height , y, of 120 inches; and thickness, T, of 
3 inches, we shall call for ten nodes in the x-direction, five tubes two 
inches wide in the y-direction, and two structure lumps between tubes, 
resulting in the nodal breakdown shown in Figure 4-6. Further, we shall 
specify four incident heat zones in both the x- and y-directions , from 
which SUBANS will set up the zone configuration shown in Figure 4-7. 

Tube and Lump Type Numbering 

For serpentine flow panels, only one tube is created; for 
parallel flow panels , n , tubes are set up, the lowest numbered tube 
being nearest the top ’Of the panel. The single serpentine tube or the 
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FIGURE U-U RECTANGULAR PANEL RADIATOR MODEL 







fluid/tube lump number 



FIGURE U-6 RECTANGULAR PANEL EXAMPLE PROBLEM - NODAL BREAKDOWN 
AND DIMENSIONS 



30" Cry p.) 



FIGUEE i|-7 RECTMGULAR PANEL EXAMPLE PROBLEM - INCIDENT HEAT 
ZONE BREAKDOWN AND DIMENSIONS 









first parallel flow tube in each panel will be numbered one greater than 
the last tube on the previous panel. For the first panel in the system 
the first tube will be Number 13. 


For eelch panel being subdivided only one lump type is required 
for fluid and tube nodes. The latter calls for conduction, to the structure 
lumps above and below, and an input code is provided whereby the user may 
call for inclusion of longitudinal conduction between tube lumps. Four 
structure lump types are created for each panel, the lowest type number 
being one greater than the last previous type for the system. The first 
of these types conducts to the structure lump immediately below (negative 
y-direction) and to the structure lump immediately to the right (positive 
x-direction) ; from the example problem, nodes of this type are 11-19, 

31-39, 51-59, and 71-79* The second type conducts only to the structiu'e 
lump below it (for our example, lumps 20, UO, 6o, and 80); and the third 
type, only to the Ivunp to the right (nodes 1-9, 21-29, 4l-49, 6I-69, 8I-89, and 
91-99). Finally, the fourth type conducts to no lumps, examples being 
lumps 10, 30, 50, 70, 90, and 100. The fluid and tube type numbers, and 
the lowest of the structure lump type n\mribers will each be one greater than 
the last previous type for the system. 

Lump Numbering 


. SUBANS will divide the panel described above into nodes numbered 
as shown in Figure I4-6. From this it may be seen that the total number of 
fluid/tube lumps created per panel is 

’^fluid/tube ~ \ \ube ^ 5 = 50] 

where n is the number of nodes in the x-direction and n^ . is the number 
of tubes in the y-direction. Similarly, the totail number of . structure lumps 

is . ■ . 


n . - . = .n n n. ' [=10 x 2x5 = 100] 

structure x.y tube 

where n is the. number of structure lumpS: between tubes . 


Lump Size 


In the x-direction, the size of all lumps will be 
size^ = x/n^ [=120 inches/10 = 12 inches] 


where X is the panel length, 
by the input width, size 

^tube 

dimensions is 


In the y-direction, tube lumps are dimensioned 
[ = 2 inches], while the. size of structure 


size 


= (y - 


n 


structure 


tube 


size 


)/n 


tube 


y 



120 inches - 5(2 inches) 
10 


= 11 inches ] 


where y is t.he prinel heif;ht. The thickness of all tube and structure lumps 
is set to the same input value of thickness. The cross-sectional area of 
the fluid lumps is a direct input item, and as such is not calculated. 


Upstream and Downstream Lump Numbers 

For serpentine flow radiators , flow will be assumed to go from 
left to right in the first tube, to go from right to left in the second, tube; 
and thence, to alternate in direction from one tube to the next. Accordingly, 
fluid upstream and tube downstream lump numbers are typically set up as shown 
in Table 4-1. 

Incident Heat 

To provide incident heat information, zones will_be set up 
across the panel, with an incident heat curve number being specified for 
each zone. The number* of zones in the x- and y-directions are also input; 
and to maintain generality, the incident heat zone boundaries need not 
coincide with node boundaries. When portions of more, than one zone fall 
on a node, a single. incident heat curve is selected according to the method 
described below. 

The size of the incident heat zones . iii the x-direction is 

size^ = x/n [ = 120 inches/4 = 30 inches] 

zone zone 

where n is the number of zones in the x-direction. Similarly, in the 

zone 

y-direction 


size 


= y/n , 


[ = 120 inches/4 = 30 inches] 


zone 


zone 


These sizes must be as large or larger than the maximum lump dimension 
in their respective direction. In the case where all of a tube or structure 


1 lim-rv 1 -1 O ri T,T-J + O v% A 


zone. When most of the area of a lump lies within a single. zone boundary, 
it too takes the curve number of that zone. When the area of a liamp is 
split equally between two zones, the following order of precedence is 
followed: if the two zones are side by side, the curve number for the 
zone to the left is assigned; if the zones are above and below each other, 
the curve for the uppermost zone is used. For a lump whose area is 
equally divided among four zones, the zone up and to the left prevails. 

As an example of the use of these r\iles , we shall select several lumps 
as follows. Structiire lump 1, being entirely in zone 1, is assigned the 
incident heat curve number for zone 1. Seventy-two square inches of 
structure lump 21 are in zone 1 and sixty square inches are in zone 5; since 
zone 1 predominates,’ its curve number is chosen. Structure lump 3 has 
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TABLE4 -1 


EXAMPLE (SERPENTINE) PROBLEM 
TYPICAL UPSTREAM AND DOWNSTREAM LUMP NUMBERS 


Fluid/ Tube Lump Upstream Downstream 


1 


0 

2 

2 


1 . 

3 

9 


8 

10 

10 


9 

20 

11 


12 

21 

12 


13 

11 

19 


20 

18 

20 


10 

19 

21 


11 

22 

22 


21 

23 

29 


28 

30 

30 


29 

40 

50 


49 

0 

For parallel flow radiators, 

flow is assumed to 

flow from left to right in 

all tubes; accordingly, 

fluid upstream and tube downstream lump number; 

are typically set up as 

shown 

in Table 4 -2. 




table 4-2 


EXAMPLE (PARALLEL) PROBLEM 

TYPICAL UPSTREAM AND DOWNSTREAM LUMP NUMBERS 

Fluid/ Tube Lump 

Upstream 

Downstream 

1 


0 

2 ■’ 

Z 


1 

■ 3 

9 


8- 

10 

10 


9 

0 

11 


0 

12 

12 


■ 11 

13 

19 


18 

20 

20 


19 

■ b 

■ , , " -21 ■■ . 


0 

22 

22 


21 

23 

29 


28 

30 

30 


29 

0 

41 


0 

■ ■ - 42 

42 


41 

43 

49 


* 48 ; 

50 

50 


49 

0 



sixty-six sqiiare inches in zone 1 and sixty-six square inches in zone 2; 
therefore, the curve for the zone to the left — zone 1 — is assigned. 
Similarly, the area of tube lump 25 is exactly halved between zone 6 and 
10; the former being uppermost, its curve is used. Finally, tube lump 23 
has equally portions of its area in zones 5, 6, 9, and 10; the rule is to 
choose up and/or to the left, hence zone 5 is selected. 


U.3.1.5 Circular Model Subdivision 

The performance of SUBRAD will be illustrated by showing what 
it will do with a sample data set. Assming a panel radius, R, of 20 inches 
(Figure H-8); and thickness, T, of one inch, we shall call for four tubes 
two inches wide in the radial direction and two structure lumps between tubes. 
Furthermore we shall set the limits of tube lump length at a maximum of 5,0 
and a minimum ofl.9 inches, resulting in the nodal breakdown shown in 
Figure U-9* We shall specify four incident heat zones, from which . SUBRAD 
will set up the zone configuration shown in Figure 4-10. 

Tube and Lximp Type Numbering 

All tubes on a circular panel are assumed to flow in parallel 
with the innermost tube being the lowest numbered. For each panel, being 
subdivided, the number of fluid and lump types will be equal to the number 
of tubes, with types being the same for a given tube. Tube lumps call 
for conduction to the structure lumps above and below, and an input code 
is. provided whereby the user may call for inclusion of longitudinal con- 
duction between tube limips. The number of structure lump types created for 
each panel will be equal to the number of tubes times the number of structure 
lumps between tubes. Each type conducts to the adjoining structure lump 
in the negative © direction (See Figure U-8) and to the adjoining lump 
in the outward radial direction. Structure lumps adjacent to and inside the 
radius of a tube l\amp will therefore conduct to only one lump. From the 
example problem; type 1 conducts to only one lump (ixrnips 1-8), type 2 
conducts to two lumps (9-16), type 3 conducts to only one Iximp (17-2U), etc. 

Lump Numbering 


SUBRAD will divide the panel described above into nodes nxjmbered 
as shown in Figure 4-9* From this It may be seen that the total number of 
fluid/ tube lumps created per panel is; 

n . ^ NT . : ■ - 

fluid/tube NL(t) [ -. = 8 + I6 + I6 + 32 = 72] 

t = 1 • 


where NL(t) is equal to the' number of lumps for tube number t. The values 
of NL(t) are either directly input into the program or calculated from a 
maximum node length specified by the user. The total mmiber of structure 
lumps can be calculated as follows : 


LDIV(i) = 


NTL(i) 

NTL(l) 



ii _ P - 16 
8 ■ ~ T 



m-1 


structure 


= NTL(l) I (NS+1) + .(NS-l)[LDIV(ia)l + 2 NS[LDIV(J )] | 

^= 1+1 ■' 
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FIGURE k-Q. CIRCULAR PANEL EXAMPLE. PROBLEM - COORDINATE 
SYSTEM AND OVERALL DIMENSIONS 


a 
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FIGUHE U-9 CIRCULAR PANEL EXAJ4PLE PROBLEM- 
NODAL BREAKDOWN AND DIMENSIONS 








INCIDENT HEAT ZONES 


= 8 I (2+1(1) + (2-l)(lt) + 2 (2+2) I =120 

where: 

i = tube number 

NTL(i') = number of tube lumps for tube number i 

NTL(l) = number of tube lumps for first tube 

NS = number of structure lumps between tubes 

ra .= tube number for last tube 

Lvimp Size 

The size and number of nodes created for a circular pianel are 
governed directly by the number of tube lumps in each tube. With reference 
to this latter nimiber there are some basic groundrules for panel subdivision 
that SUBRAD follows. These include: 

(1) The number of tube lumps per tube can be specified by the 
user or calculated by SUBRAD. 

(2) User specification; starting with n lumps for the first 
(inside) tube, each successive tube is required to have 
some multiple of n lumps; i.e.,ln,2n,3n, etc. 

(3) Automatic SUBRAD calculation; the user specifies a set of 
limits for tube lump length which SUBRAD will use to determine 
the number of lumps per tube, attempting to stay within the 
specified limits for each lump, while following the guideline 
as described above. If such a breakdown is not possible 
SUBRAD will override the maximum length specified by the 
user. 

(U) Structure lumps inside the radius of the first tube-coin- , 

■ Cide with the. theta angle (Figure 1+-8) of tube lumps 
in the first tube (See Figure 4-9). 

. (5) All other structure lumps coincide with the theta angle 

of the tube lump adjacent to and nearest the center of the 
panel. 

Upstream and Downstream Lump Numbers 

For circular radiator panels, flow is assumed to flow clockwise 
in all tubes; accordingly, fluid upstream and tube downstream lump numbers 
are typically set up as shown in Table 4-3. 

Incident Heat 


To provide incident heat information, zones will be set up 
across the panel i with an incident heat curve number being specified for 
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TABLE It- 3 

.CIRCUI.AR PANEL -EXAMPLE PROBLEM 
TYPICAL UPSTREAM AND . DOWNSTREAM LUMP NUMBERS. 


FLUID/TUBE LUMP 

1 

2 

7 

8 
9 

10 

23 

2k 

25 

26 

55 

56 

57 

58 

87 

88 


UPSTREAM 

0 

1 

6 

7 

b 

9 

22 

23 

0 

25 

5U 

55 

0 

57 

86 

87 


DOWNSTREAM 


3 

8 

0 : 

10 

11 

2lt 

0 

26 

27 

56 . 

0 

58 

59 
88 

0 
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each zone. The number of zones on the panel is input; and to maintain 
generality, the incident heat zone boundaries need not coincide with node 
boundaries. When portions of more than one zone fall on a node, the 
zone covering the greatest portion of the area will be assigned to the 
node. When the area of a lump is split equally between two zones the 
order in which the zones are specified determines which zone is assigned 
to the lump. The node will always fall in the zone with the highest order 
but will never fall in the last zone. To clarify this point, consider a 
panel with 6 zones, the first zone being of order 1 and the last of order 6. 
If a lump is split between zones 3 and U, U will prevail; for a split 
between 1 and 6 , 1 will prevail; for a split between 5 and 6 , 5 will prevail. 

For circular panels incident heat zones are determined by 
boundary lines, input by the user, by specifying the x and y-axis intercepts 
of the line. If a boundary line goes through the origin of the panel co- 
ordinates, a theta angle (see Figure H-8) must be specified. Zones may 
be specified by any number, of boundary, lines , .yielding: the., ability to des- 
cribe almost, any size,, shape or form of zone by proper use of these lines . 
Boundary lines are input in sets of one or two along with an incident 

heat zone number. 

SUBRAD’s interpretation of these boundary lines can best be 
described by the groundrules built into the subroutine as listed below: 

( 1 ) For a panel with only one incident heat zone, boundary 
lines are not required. 

( 2 ) For a zone described by only one boundary line, that 
particular zone will be the larger of the two created by' 
the boundary. If the panel is split into two zones of 
equal area, the zone which includes the lowest numbered 
node will be assigned to that particular incident heat 
ciurve number . 

(3) For a zone described by two boundary lines, the. nodes 
falling between the two boundaries will be assigned to 
the zone. If the two boundary lines intersect at some 
point, the zone described by these boundaries is determined 
as shown in Figure 4-11. 

(4) For zones created by more than one set of boundary lines, 
the following procedure is followed; SUBRAD sets up pre- 
liminary zone numbers for the first set of boundary lines; 
any successive set of boundaries will override a previously 
assigned nodal zone number specified by a. previous set 

of boundaries. 

( 5 ) Any nodes not falling within the prescribed boundary limits 
of any of the incident heat zones are assumed to lie in a 
zone with an incident heat curve number indicated for the 
last zone . 
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As an example of these rules, we will select several nodes as 
shown on Figure 4-10 of our example problem and describe how SUBRAD inter- 
prets the information supplied by the user. Three sets of boundary lines are 
necessary to describe the four incident heat zones on our example panel. 

Table 4-4 gives a description of the incident heat zone information supplied 
by the user. The first card being of order 1 describes lines AA and BB of 
Figure 4-10. All nodes falling between these lines are assigned incident 
heat curve number 10. Structure lumps 51j and 85 do not lie entirely 
between the boundaries but over half of their area does, therefore these 
lumps are assigned to zone number 10. Boundary lines BB and CC describe zone 
number 20 having an order of 2, while CC and DD describe zone number 30 of 
order 3. Even though structure lumps 101 and 117 are split evenly between 
zones 20 and 30, both lumps are assigned to zone 30 since its order of 3 is 
higher than that of zone 20 whose order is 2. The remaining lumps which do 
not fall in any of the above 3 zones are assigned to the zone specified by 
the user as the last incident heat zone. For our example the zone is number 
40, but could have been 10, 20, 30 or any number so desired by the user. 

4. 3.1.6 Two-Dimensional Model Subdivision 

The performance of SUB2D will be illustrated by showing what it 
will do with a sample data set. We shall assume for the example a panel with 
four tubes and two structure lumps between tubes. Furthermore we shall set 
the number of tube lumps per tube and let SUB2D divide the model accordingly. 

Figure 4-12 shows how the panel is subdivided into lumps. Lumps 
are created by lines originating from a central point and dividing each tube 
length into lumps of equal length. Twelve lumps were prescribed for the 
first two tubes and 24 lumps for the last two tubes. Subsequent divisions 
can be made in order to reduce the lump len^h for any tube past the first as 
in tube number three of the example 

Automatic Tube Length Option 

The number of lumps per tube can be prescribed by the program 
user, but in addition an automatic tube length option exists in SUB2D. This 
option allows the program user to specify a set of limits of tube lump length 
and SUB2D will automatically subdivide the panel accordingly. For this 
automatic division the user indicates the maximum length for a tube node in 
the first (innermost) tube whereupon the program sets this length at a value 
less than the prescribed value (XMAXl). An additional limit on node length 
(XMAX2) for the rest of the tubes is provided by the user also. XMAX2 is 
the maximum length for tube nodes past the first tube. With the divisions 
set up by the first tube^ the routine checks the value of the tube node 
length for tube number two. If this length is greater than XMAX2, then a 
subsequent division is made making twice as many lumps in tube 2 as in tube 
i. This process is continued for the entire panel. 

For panels with a large number of tubes, the number of lumps can 
become excessive if very many subsequent divisions past the first tube are 
allowed. For this reason an additional check is made on the parameter NDIVA, 
another user input., This parameter is used to limit the number of lumps per 
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TUBE . 


TUBE #3 


TUBE #2 


TUBE §1 



FiGuiffi h -12 two-dimensiq:^al panel nodal , breakdown 


panel. This number represents the maximum number of subsequent divisions 
allowable past the division in the first tube. If after a test of tube 
lump length against X^X2, another division is necessary^ the routine checks 
to see if this would exceed the limits set up by NDIVA. If so, no further 
division in the panel is allowed. 

Variable Tube Spacing Option 

SUB2D has a variable tube spacing option whereby the tube spacing 
for the first tube can be different than that of the other tubes. By setting 
TSPl (Card ANS-ID) to the desired value a unique tube spacing is set up for 
the first tube and the remaining tubes are s^ced equally. If equal tube 
spacing is desired for the entire panel TSPl should be set equal to zero or 
just left blank. 

Tube and Lump Type Numbering 

All tubes on a two-dimensional rectangular panel are assumed 
to flow in parallel with the innermost tube' be ii^ the 'lowest ’numbered . Tube 
lumps call for conduction to the structure lumps above and below, and an 
input code is provided whereby the user may call for inclusion of longitudinal 
conduction between tube lumps. Each structure lump conducts to the adjoining 
structure lump in the outward radial direction and to the adjoining structure 
lump in the positive 9 direction (See Figure 4-12). Structure lumps adjacent 
to and inside the radius of a tube lump will therefore conduct to no more 
than one lump. 

Lump type numbering can best be understood by referring to 
Figure 4-13. Tube lump types start at 0 = 0° position and increase 
numerically to 0 = 45° position whereupon the type number: decreases by one 
to 0 = 90 ° position. The process is then repeated for the rest of the 
panel. 

Structure lump types also start at 0 = 0° position and then in- 
crease numerically to 0 = 90° position. This process is then repeated for 
the rest of the panel except for the last lump. Since the last lump does 
not conduct to any other lump it is assigned a new type number. 

Lump Numbering 

SUB2D will divide the panel described above into nodes numbered 
as shown in Figure 4-i4. Lump numbers start at the 0 = 0° position and 
increase numerically around the panel. Tube lumps and structure lumps, are 
allowed to have identical numbers since the identification of each is known 
internally in the program. 

Lump Size 

The size and number of nodes created for a two-dimensional 
panel are governed directly by the number of tube lumps in each tube. With 
reference to this latter number there are some basic ground rules for panel 
subdivision that SUB2D follows. These include; 
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STRUCTURE TYPE 
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FIGURE 4-13 TWO-DIMENSIONAL- PANEL LUMP TYPE CLASSIFICATION 


















STRUCTUHE LIMP 
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PIGURE h~lk TOD-DIHENSrONAL PANEL LUMP NUMBERING 


















(1) The number of tube lumps per tube can be specified by 
the user or calculated by SUB2D, 

(2) User specification; starting with n lumps for the first 
(innermost) tube, each successive tube is required to have 
some multiple of n lumps; l.e., In, 2n, 3n, etc. 

(3) Automatic SUB2D calculation; the user specifies a set 

of limits for tube lump leh^h whereupon SUB2D will deter- 
mine the number of lumps per tube attempting to stay 
within the specified limits for each lump, while following 
the guidelines as described above. 

(4) Structure lumps inside the radius of , the first tube 
coincide with the theta angle pf tube lumps in the first 
tube. 

(5) All other structure lumps coincide with the theta angle 
of the tube lump adjacent to and nearest the origin of 
the panel. 

Upstream and Downstream Lump Numbers 

For 2-D radiator panels, flow is assumed to flow counter-clock- 
wise in all tubes; accordingly, fluid upstream and tube downstream lumps 
numbers are typically set up as shown in Table , 4-5 v 

Incident Heat 

Only one incident heat zone is allowed for the entire panel. 
4.3.2 Plot Options 

A positive integer punch in column 72 of parameter card 2 will 
cause the generation of a plot file. This tape will have ail of the fluid, 
tube, and structure lump temperatures, plus other items indicateii below. 

The plot file is generated bn tape unit I. 

The format of the plot tape is; 

Record No. 1 

Title (from title card, 12A6), 0, 0, 0, 0, 0, 0, 1, 0, 0, number 
of pressure drops, ,1, number of flow rates, number of fluid 
temperatures, number of tube temperatures, number of structure 
temperatures. 

Record No. 2 

Time, total pressure drop, tube pressure drops, heat rejection, 
flow rates, fluid temperatures, tube teniperatures, structure 
Temperatures. 



SUBRAD CARD NUMBER RAD- 
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TABLE k-h CIRCULAR PANEL EXAMPLE PROBLEM - 
INCIDENT HEAT ZONE SPECIFICATION 































WO-DIMENSIONAL PANEL EXAMPLE PROBLEM 
TYl’ICAL UPSTREAM AND DOWNSTREAM LUMP NUMBERS 
(REFE^NCE FIGURE it-B) 


Fluid/Tube Lumt 



Downstream 


1 

O' ■ 

2 

2 

■ 1 ■ ■ 

3 

7..- • 

6 

8 

11 

10 

12 

12 

11 

0 

13 

0 

14 

14 : 

13 

15 

19 

18 

20 

23 

22 

24 

24 

■ ' 23'-:' 

0 

25 

.0" 

26 

26 

25 

.27 

37 

. 36 

38 

47 

. . 46 

48 

48 

.'47 -■ 

0 

49 

■■■■■. ;,0'' 

50 

50 

49 

51 

6l 

60 

62 

71 

70 

72 

72 

71 

0 




Record No. 3 

Time , etc . as on Record 2 

Tlif! Uiol. l•eco^i him i.i tieK.'<'ti v‘) Umu to huili'Hl.r ^<iui oT uni. put. 
U.3.3 Checkout Print 

It is possible to obtain a detailed print at each iteration (this 
is usually used in data checkout* and so is called checkout printing). 

This can be obtained by putting a "l" punch in column 60 of parameter 
ca,rd 2 . The checkout print provides a direct, relatively simple check 
of the data input consistency, and indicates, the minimum calculation time 
increments calculated for each lump. A survey, of these time increments 
permits choice of a minimvim time increment which .will acourately character- 
ize the transient and override only these lumps which can realistically 
be considered at steady state. 

U . 3 . 4 Restart 

R eq,uest.ed. Pump...f or, Restarting 

Any problem can be dumped and restarted at a later time. This 
is achieved by punching a "l" in column 58 on parameter card 2 . This 
option is useful in data checkout in that a problem can. be submitted for 
a short transient time, and, after examination of the results, restarted 
for a longer transient time. . The computer req.uest card must specify that 
the, output tape is expected, and the proper set-up card must.be included 
in the deck. 

Restart Procedure 

The procedure for restarting a problem which has beendumped 

is; . 

■ Cl) Fill .'out the computer request , card as . in. an initial run, 
except specify, the previously dumped tape as an input 
tape on tape unit J. 

(2) .Submit only the first two of the data cards (that is, 
parameter cards 1. and ?} with a "l” punch ill col’jmn 62 
on Card 2 to indicate that data is to be read from a 
restart' tape. 

^• 3.5 . EDIT 

The. large number of data cards required for problems run on 
this routine presents three problems; (l) increased probability of 
operator and/or card reader error, (2) increased probability of a card 
reader jam, and ( 3 ) . significant extra time required to read in data from 
the card reader when problem (2) occurs. For these reasons a routine 
was developed for; reading input . data from tape with the capability for 
modifying the data on read-in. 

The EDIT routine is called by parameter INDATA input in • 
columns 6 j and 68 on parameter card 2. Possible inputs are: 

Cl) INDATA = 0, All data is supplied on cards. 
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(2) INDATA = 1, All data is supplied on cards and the card 
images are written. on tape on Unit B. 

(Should be specified as an output tape on Job 
card) . 

. (3)., INDATA = 2,. Use data input on tape on unit ,D with desired 

changes on cards, to write a new data tape on 
unit B. (D is input tape and B is output tape)* 

(It) INDATA = 3) Use the data, read in from unit B without change. 
(B is input tape). 

(5) INDATA = Uj List, and use data read in from unit B without 

. • change. (B is input tape). 

(6) INDATA - 3 , /^nch and use data read in from unit ,B. without 

. . change. (B is input tape). 

For INDATA =. 3 » and .5, parameter cards 1 and 2 are read in 

from cards . - 


When INDATA = 2, the deck set-up consists of pa,rameter cards 1 
and 2, the EDIT control cards (described below),- and the new data cards 
(with. the same format as the cards being replaced). 


2, are; 

The EDIT 

control cards, used 

only When INDATA has a value of 

Columns 

Format 

Nomenclature 

Description 

1-5 


■ . K3 : ■ 

Card number K 3 will be either 
.removed or pushed down depending 
ubon, the contents of the next two 
fields (Kl4 and. K5). If K3 is 
greater than the number of cards 
on the tape plus 10000, no edit, 
will be performed. 

6-10 

15 

., k4 

kU must be equal to or greater 


than K3 or must be left blank. 

If ^ K3, cards numbered K3 

to inclusive will not be 

transferred to the hew data tape 
on. unit B. If left blank no cards 
will be deleted. 

11-15 , 15 K5 If not blank or zero, K5 cards, 

which must immediately follow this 
.edit control card will be inserted 
ahead of card K3 or in place of 
cards K3 through kU. 

Giving K3 the value of 99999 will cause the continuing transfer of card 
images from Unit D '.to Unit B until a card image , containing bbbl3b ,(b = 
blank) in columns 1-6 is found. The Unit D tape is .not altered in any way 
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should there be errors in the edit, deck which cause fatal errors in the LTV 
program. It is the responsibility of the user to maintain extra copies of 
the data tape and/or an up-to-date card deck. 

If automatic nodal subdivision is used in a run. Unit C must be 
substituted for. Unit B in the above discussion of the EDIT routine. The data 
on cards or Unit C is combined with output on Unit E from. automatic nodal 
subdivision and resulting data is placed on Unit B. 

U.3.6 Specified Temperature Nodes 

If the mass-specific heat product for any node is zero, the 
temperature of that node will remain constant at the initial temperature 
throughout the dxiratioh of the problem. 

A temperature history may be imposed on either structure or 
tube nodes for the duration of the problem. If a node has an imposed 
temperature variation, the prescribed temperature curve must be input 
for the node. 
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4.4 . • User’s Manuals 

4 .4.1 Transient Performance Routine (TPR) 

The user’s instructions for the Transient Performance Routine 
are presented in this section. . 

Parameter Cards 


Columns 

Format 

Nomenclature 

Description 

Parameter 

Card 1 



1-72 

12A6 

: HEADER ' , 

Any 72 alphanumeric characters 
to be used for page heading. 

Parameter 

Card 2 



1-10 

FIO .5 

TIME 

Transient start time, hr. 

11-20 

FIO .5 . 

TAU ' 

Transient stop time, hr. 

21-30 . 

FIO .5 

TINCMN . : 

Time increment, hr. 

31-35 

F5.0 

, DELTAU 

Print interval, hr. If 
DELTAU < TINCMN, print 
interval = TINCMN. , 

36-4o . 

F5.0 

RTIME 

Computer time requested, minutes. 

4i-45 

F5.0 

THETA •’ 

"Convergence factor-. Routine sets 
to .9 if left blank. 

46-50 

F5.0 

DPTOL ' , 

Parallel flow, path pressure drop 
balance tolerance (decimal 
fraction) . , 

51-55 

F5.0 

SSTEST 

Steady state tolerance, °F. 

56-58 

12 

IDUMP 

= 0, No dump tape to, be written. 
=1, Dump data on UNIT I when 
either TAU or RTIME is 
exceeded. 

59-60 

12 

NCKOUT 

= 0, No checkout print. 

= 1 , Checkout, print . ' 

61-62 

' 12 ■ 

isTart 

= 0, New, data follows. 

=1, Read data from restart tape. 

63-64 

12 

NVLVRS 

Number, of valves having parameter 


to be changed on restart. 
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Columns 

Format 

Nomenclature 

65-66 

12 

NORAD 

6t^68 

12 

INDATA 





■ 

69-72 

I4 

NPLOT.; 



Parameter 

Card 3 


1-5 

15 

NTUBE. 

6-10 

15 

NFLT. 

11-15 

15 

NMLT • 

16-20 

15 

NT 

21-25 

15 

■ NTFL 

26-30 

.15 

NTML 

31-35 

15 

NSL „■ 

36-Uo 

15 

NPSl 

41-1+5 

15 

NPS2 

Parameter 

Card 4 

(Leave out for NORAD 

1-5 

15 

. NTYPE(l) 

6-10 

15 

. NTYPE(2) 


Description 

Number of radiator panels to be 

divided into nodes . 

= 0, All data supplied on cards. 

=. 1 , Write data cards on UNIT B 
(UNIT C if NORAD > O). 

= 2, Use data on UNIT D plus 

edit cards to create data 
‘ on UNIT B (UNIT C If NORAD >0) 

= 3, Use data on UNIT B (UNIT C 
if NORAD >0). 

= U, List and use data on UNIT B 
(UNIT C if NORAD > O). 

= 5 j Punch on cards and use data 
UNIT.B- (UNIT e if NORAD > 0). 

^ 0, Plot file output will occur 
every NPLOT iteration on 
UNIT I. ' 

= 0, No plot file will be generated 
■ on UNIT I. 


Number of tubes . • 

Number of fluid lump types . 

. Number of tube lump types . 

Number of structure l\mip types. 

Number of fluid Itamps . 

Number of tube lumps . 

: Number of structure limps . 

Number of ■ panels on Side 1. 

Number of panels on Side'. 2. 

= 0 ) ■ . . 

= 1, Rectangular panel. 

= 2, Circular panel. 

® 3> 2-D Rectangular panel. 
Panel type for second panel on 
Side 1 . or panel on side 2 
if only 1 on side 1. 



Columns 

Format 

Nomenclature 

Description 

,11-15 

15 

NTYPE(N) 

Panel type 1'or.iaut patiei (jti 

I- 

16-20 

15 

’ .NTYPE(N+ 1 ) ■ 

Panel type for first panel bn Side 
2 -' 

21-25 . . 

15 . 

ntypeCm) . 

Panel type for last panel on Side 

, 2 .. • ' 

Insert ANS data dards 

here if .NORAD > 0 

(Pa.rameter Card 2, Columns 65-66). 

Parameter 

Card 5 



1-5 

F5.0 

•HIl . .. 

Entry length heat transfer co- 
efficient factor. Recommended 
value for circular tubes is 0 . 575 - 

6-10 

F5..0 

HI2- 

■Fully developed heat transfer 
coefficient factor. Recommended 
value for circular tube is 1 . 0 . 

11-15 

F5 .0 

ANEW ... 

Fraction of new flow rate to use 
when averaging flow rates. 
Recommended value is 0.5- 

16-20 . 

15 

. NGYCLE . 

Time dependent curve option. 
- 0, Non-cyclic. 

= 1 ,. Cyclic . 

20-25 

15 . 

KLSDLP 

Closed loop option 
= 0, Inlet temperature to tube 1 
supplied on curve (Parameter 
Card 8, Cols. 31-35) 

= 1, Outlet temperature of tube 
8 used as inlet temperature 
to tube 1. (Parameter Card 8, 
Cols. 31-35 may be left 
blank ) . 

Parameter 

Card € 



.1-5 

.15 

NTl 

*Nvimber of tubes in panel 1. 

6-10 

15 , 

NT2 

Humber of tubes • in panel 2. 

11-15 

15 

NTS 

Number of tubes in panel 3. 

Continue 

listing the 

ntmiber of tubes in 

each panel in five column fields. 


* Set number of tubes to a minus quantity for: circular arid 
two-dimensional panels . 
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Parameter Card 7 


Coliunn 

Format 

Nome nc la t lire 

Description 

1-5 


LMPINl 

t) 

Inlet lump for side 1 of the 
regenerator. 

6-10 . 

15 

LMPOTl 

Outlet lump for side 1 of the 
regenerator. 

11-15 

15 

LMPIN2 

Inlet Jump for side 2 of the 
regenerator. 

16-20 

15 

LMP0T2 

Outlet liomp for side 2 of the 
"^regenerator.” ^ — 

21-30 

F10.5 

UA 

Product of the overall heat 
transfer coefficient and the area 
of heat transfer (UA) for the 
regenerator, BTU/hr-°F. 

Parameter 

Card 8 



1-5 

15 

. NDENC 

Fluid, density curve number. 

6-10 

15^ 

NCONC , 

Fluid conductivity curve number. 

11-15 

15 

NSHC 

Fluid specific heat curve niimber. 

16-20 

15 

NVISC- 

Fluid viscosity curve n\omber. 

21-25 

15 

NFFC 

Fiuld friction factor curve number 
for turbulent flow . 

26-30 

15 

NFLOW 

Inlet flow rate curve number. 

31-35 

15 

ntemp 

Inlet temperature curve number. 


Ten cards are required to specify the valve parameters. Valve 
parameters for the radiator bypass valve (valve number 1) and the regeneratoi' 
bypass (valve number 2) are specified on two pairs oV Carti A and Cnvd lU 
Valve parameters for the proportioning valve (valve number 3) are specified 
on Cards C and D. Valve parameters for the two polynomial or rate limit 
stagnation valves (valve n\xmber . 4 and valve, number . 5) are specified on two 
pairs of Card A and Card B.. Valves do not need to be in. numerical order. 


Bypass and Stagnation Valves 

Format Nomenclature Description 

13 NVN Valve niomber. 

12 NOP Operation mode, * 

. = 0, valve is operating. 

=1, valve is locked in initial 
. position. 

=2, Valve is not considered; i.e. 
flow is distributed among 
all tubes on panel according 
to pressure drops (for valves 
4 and 5 only) . 


6-10 

15 

NSLN . 

Sensor lump number. 

11-20 

FIO . 5 

FRMIN 

Minimum fraction 
non-bypass tube . 

allowed through 

21-30 

F10.5 

FRMAN 

Maximum fraction 
non-bypass tube . - 

allowed through 

31-40 

FIO. 5 

FRIN 

Initial fraction 
non-bypas s tub e . 

allowed through 


Card A - 
Columns 
1-3 
4-5 



Columns Format Nomenclature Description 

Card B 


1-5 

15 

NTYPE 

= -1, Bypass valve for solar 
a*bsorber . 

=1, Rate limit bypass. 

= 2, Polynomial bypass . 

If NTYPE 

II 

1 

H 



6-15 

FIO . 5 

SETPT 

Set point temperature, °F. 

16-25 

F10.5 

DBAND 

Dead band, °F. 

26-35 

FIO . 5 

REACT 

Rate factor. Units are fraction 
-= bypass per -second -per ^F. 

36-ii5 

FIO. 5 

RLIM 

Rate limit. Units are fraction 
bypass per second. 

If NTYPE 

= 2,-1 



6-15 

FIO. 5 

> 

o 

Coefficients in fraction bypass . 

16-25 

FIO . 5 

^1 

X = Aq+A^T+A^T^+A^T^+A^T** 

26-35 

FIO. 5 

OJ 

< 

where T is sensor lump tempera- 
ture (°f). 

36-U5 . 

FIO . 5 



46-55 

FIO. 5 



Card C Proportioning!: Valve 



1-3 

13 

NVN 

Valve numher. 

1, c 

TO 
-1. £- 

■MOTD 

^ ^ VS ^ 

=0, Valve is operating. 

= 1, All the radiator flow is 
directed through side 2. 

6-10 



Blank 

11-20 

FIO. 5 

POSMIN . 

Minimum allowable linear valve 
position from left, milli-inches 

21-30 

FIO. 5 

POSMAX 

Maximum allowable linear valve 


position from, left , milli-inches 
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Columns 

Fonnat 

■ Nomenclature 

Description 

31-Uo 

FIO .5 

POSIN . 

Initial linear valve position from 
left, milli-inches (mils). 

Card D - 

Proportioning Valve 


1-10 

. FIO . 5 . ■ 

. . FULOPN . 

Maximum linear valve position 
from left, mils. 

11-20 

FIO .5 

VLVGAN 

Valve gain, miili-inches/°F. 

21-30 , 

FIO. 5 

PPARA 

Panel parameter. Enter 2. 

31-i+O 

FIO. 5 

. GFACT . 

Geometry factor. 2 

GFACT .= (2.92)(1D”°)(D )* 
where D is valve orifice diameter 
in mils. Units are 




(is.) 2 (-^)^ ( ■ r-) 

'hr/, 'rail'^^ ''1000 psi' 

1+1^50 

Flo. 5 

VLVTOL 

- N\ill position tolerance, mils. 

51-60 

FIO. 5 

. TCON 

Time constant, seconds. 


Fluid Data 

Cards 


Card 1 




1-12 



Blank 

13-36 

1*A6 

, ALPHA . 

FLUID b TYPE b DjITA bbbbbbbbb 
where b denotes a blank. 

Card 2 




1-25 



Blank 

26-30 

15 

NKPDC . 

Dynamic head loss, K. 

Factor, curve ;n\imber (leave blank 
if not needed) . . 

31-ltO 



Blank 

Ul-50 , 

FIO. 5 

. FLL 

Fluid lump length, in. 

51-60 ^ 

FIO. 5 

CSA 

Cross-sectional area, sq.in. 

61-67 

FT -5 

WP 

Wetted perimeter, in. 

68-72 

F5.it 

FRE 

Factor for computing friction 


factor as a function of Reynold's 
niunber. Routine sets to 1.0 if 
left blank. 

* Value for Apollo Block II Environmental Control System Valve. 


Repeat Card 2 for every fluid type. 
Card 3 

Columns Format . Nomenclature 

1-12 


13-36 

.4a6 ; 

ALPHA 

Card 4 



.1-5 

15 

LN 

6-10 

15 

NLU 

11-15 

15 

NTB 

16-20 

15 

NTYP 

21-30 

FIO . 5 

TI 

Repeat Card 4 for every fluid lump. 
NTFL and must be in niomerical order 


Tube Data 

Cards 

Card 1 



1-12 



13-36 

ltA6 

ALPHA 

Card 2 



1-5 

F5.2 

DEN 

6-10 

15 

NCONC 

11-15: 

15. 

NSHC 

16-20 

15 

NABC 

21-25 

15 

NEMC 

26-30 

15 

NTCT 

31-35 

15 

. NFCT , 


Description 

Blank 

FLUID L LUMP b DATA bbbbbbbbb 
where b denotes a blank. 

Lump number. 

Lump upstream. NLU=0 for first 
lump 3 ji_eyery, ^i^e. 

Tube . number . 

Type number. 

Initial temperature, °F. 

The l\amps must be numbered 1 through 


Blank 

TUBE b TYPE b DATA bbbbbbbbb 
where b denotes a blank. 

Density, Ib/ff^. 

Conductivity curve number. 

Specific heat curve nvunber. 

Absorptivity curve number. 

Egiissivity curve number. 

Number of tube lumps conducted 
"to". 

Number, of structure lumps conducted 
"to". 

96 . . 



Columns 

Format 

Nomenclature 

Description 

36-Uo 

15 

LCC 

= 0 , Lumps of -this type have 
longitudinal conduction. 

= 1, Lumps of, this type do not 

have longitudinal conduction. 

Ul-50 

FIO .5 

XI . 

Dimensions, in. X 1 *X 2 *X 3 = volume 

51-60 

FIO .5 

X2 

Xl /2 = longitudinal conduction distance. 

61-70 

Card 3 

FIO .5 

X3 

X 2 ;X 3 = area for longitudinal conduction 

1-10 

FIO .5 

ART 

Area for heat : transfer to enclosed 
fluid Imp, sq. in. 

11t-20 

FIO .5 

.ae: 

. Surface area for external radia- 
tion,, sq. in. 

21-30 

FlO.5 . 

RAC 

Factor for dividing conduction 


distances and dimensions. 

Routine sets to 1.0 if left blank. 


Use when conduction dimensions 
do not fit into the five column 
fields. 


31-35 

F 5.5 

Y 1 

Conduction distance of this 
type of ’*frbm" lump to first ”to*^ 
lump . 

36-uo 

F 5 . 5 

Y 2 

Conduction distance for first ”to” 
Imp listed for this type of 
’’from". Imp. 

kl-U 5 

F 5.5 

B 

B X D is area for conduction for 
first "to" Imp for this type of Imp 

U6-50 

F 5.5 

T) ' 


51-55 

F 5.5 

Yl. 

Conduction data for second "to" 

56-60 

F 5.5 

Y 2 

Imp. . 

61-65 

F 5.5 

B 


66-70 

F 5.5 

D . . 


If NTCT 

+ NFCT > 2, 

follow with Card 4 


Card 4 




1-10 



Blank 

11-15 

F 5.5 

Yl 

Conduction data for third "to" 

16-20 

F 5.5 

. Y 2 

Imp . 

21-25 

F 5.5 

B 


26-30 

F 5.5 

. ■ D 
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Columns 

Format 

Nomenclature 

Description 

31-35 

F5.5 

Y1 

Conduction data for fourth "to" 

36-40. 

F5.5 

Y2 

lump. 

4i-45 

. F5.5 

. ■ B ■ 


46-50 

F5.5 

. D 


51-55 

F5.5 

1 

Y1 

Conduction data for fifth "to" 

56-60 

F5.5 

. • Y2 

lump . 

61-65 

F5.5 

B 


66-70 

F5.5 

D 



Repeat Card k if NTCT + NFCT>5« If NTCT > 0, the data for tube lump 
"to" tube lump must.be given before the data for tube lump "to" structure 
lump . , 

Repeat Cards 2 and 3 (followed by Card 4 if needed) for every tube l\mip 
■ type. - -- - •. — ■ ■ ; - ■ - - - 

Card 5 


1-12 



Blank 

13-36 

4a6 

ALPHA' 

TUBEb LUMP b DATA bbbbbbbbb 
where b denotes blank. 

37-80 



Blank 

Card 6 

( Tube Lump Cards ), . .( One for ,e%ch 
numbered 1 through NTML and must 

tube.,l\mip. The lumps must be 
be entered in numerical order) 

1-5 

15 

■ LN 

Lump number. ' 

6-10 

15 

NDL 

Lump number of tube lump down- 
stream. 

= 0, for last lump in each tube 

11-15 



Blank 

16-20 

15 

NTYPE 

Type number of lump 

21-30 

FlO-5 

TI 

Initial temperature, °F. 

31-35 

15 , ■ 

NQIC 

Incident heat curve number. 

36-40 

15 

NTWC 

Prescribed temperature curve 
number (may be left blank) 

41-45 

15 " 

NTLl 

First lump conducted "to". 

46-50 

■ 15 

NTL2 

Second lump conducted "to". 
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Columns 

Format 

■Nomenclature 

Description 

51-55 

15 , 

NTL 3 

Etc . 


56-60 

15 

• NTLlt 

The order in which the lumps 
conducted "to" are listed will 

61-65 

15 

NTL 5 . ■ 

depend upon the order in which 
the conduction data was given on 

66-70 

.15 

. NTL6 

the type cards . The tube liunps 
conducted "to" (if any) must be 


listed before the structure lumps 
conducted "to" (if any). 


Card 7 (Continuation of list of lumps conducted "to". If the nvunber of 

liunps conducted , "to" ■ is greater than 6 j follow Card 6 with Card ?)• 


1-5 

15 . 

NTL 

Next lump conducted "to". 

6-10 

•' 15 , 

,ntl' 

Next lump conducted "to". . 

Etc . to 




66-70 

15 

NTL .. 


Repeat Card 7 as needed to list all lumps 


Repeat. Card 6 followed by Card; 7 (if required) for every lump. The lumps 
must be given in increasing rivimerical order. 


Structure. Data 

Cards 


Card 1 




1-12 



Bla;nk 

13-36 

U6 

■ ALPHA ^ 

STRUCTURE b ' TYPE b DATE bbbbbbbbb 
where b denotes a blank 

37-80 



Blank 

Card 2 




1-5 

F 5.2 

DEN 

Density of fin material, Ib/ft^. 

6-10 

15; ■ ' ‘ 

NCONC 

Conductivity curve number. 

11-15 

"15 .' 

.. NSHC 

Specific heat curve nvimber. 

16-20 

15 

NABC 

Absorptivity cvirve number. 

21 - 25 ; 

15 ' 

NEMC:. ' ; . ' ' 

Emissivity curve number.- 


9 ? 



Columns 


Format 


Nomenclature 


Description 


26-30 

• 15 

NFCT 

31-1+0 

FIO.5 

XI 

1+1-50 

F10.5 

X2, 

51-60 

FIO.5 

X3 

If NFCT 

> 0, enter Card 3. 


Card 3 



1-10 

FIO.5 

.FAC 


~ ■ 



11-15 

F5.5 

Yl 

16-20 

F5.5 

Y2 

21-25 

F5.5 . 

B 

26-30 

F5.5 

D 

31-35 

F5.5 

Yl 

36-1+0 

F5.5 

Y2 

1+1-45 

F5.5 

B 

46-50 

F5.5 

D 

51-55 

F5.5 

Yl , 

56-60 

F5.5 

Y2 . 

61-65 

F5.5 

B 

66-70 

F5.5 

D 

■ If NFCT 

> 3, repeat Card 

3 as needed 

cards which are a repeat 

of Card 3. 

Repeat 

Card 2 (followed by Card. 3 if 

Card' 4 

1-12 



13^36 

4a6 

Alpha 


37-80 


Number of fin lumps conducted 
"to" by lumps of this type. 

Dimensions of lump, inches. 
X1-X2-X3 = Volume. 

X1*X2 = External radiation area. 


Factor for dividing conduction 
distances and dimensions. Routine 
sets to 1 if not given, use when 
numbers are too small to fit five 
column field. 


Data for. conduction to first lump 
conducted "to" by this type lump. 


Data for second lump conducted 
"to" by this type lump. 


•Data for third lujtip conducted' 
"to" by this type lump. 


■ FAC should be omitted' oh all 


needed) for every structure lump type. 


Blank 


STRUCTURE b LUMP b DATA bbbbbbbbb 
where b denotes a blank 

Blank 


Card 5 (Structure Lmp Cards) (One for each structure lump. The lumps 
must be numbered 1 through NSL and must be entered , in numerical 
order) 
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Columns . 

Format 

Nomenclature 

Description 

1-5 

15 . 

LN 

Lump number*. 

6-10 

15 . 

LTYPE 

Type number of lump 


11-20 

Flo . 5 

TL 

Initial temperature of lump, 

21-25 

15 

NQIC. 

Incident heat curve number. 

26-30 

15 ; 

NTWC 

Prescribed temperature curve 
number (may be left blai^). 

31-35 

15 

NTLl 

First liamp conducted "to". 

36-40 

' ■ 15 

NTL2 

Second lump conducted "to". 

41-45 

15 

NTL3 

Etc . 

46-50 

15 

. NTL4 

The ;order in which the lumps 

51-55 

15 

NTL5 

conducted "to" are listed will 

56-60 

15 . 

NTL 6 

depend upon the order in which 

61-65 

15 

. NTL7 

the conduction data was given 

66-70 

15 

NTL 8 

on the type cards . 


Card 6 (Continuation, of list of lumps conducted ”to". If the number of lumps 
conducted ”to” is greater than 8 , follow Card 5 with Card 6 ). 


1-5 

15 

NTL9 

Ninth lump conducted "to". 

6-10 

15 

NTLIO 

Tenth lump conducted "to". 

11-15 

15 

. NTLll , 

Eleventh lump conducted "to" 

Etc . to 




66-70 

15 

NTL22 

Twenty-second lump conducted 


Repeat Card 6 if mmber of liomps conducted "to" is greater than 22, 

Repeat Card 5 followed by Card 6 (if required) for every lump. The lumps 
must be entered in increasing numerical order. 

Curve Data Cards 

Card 1 

1-12 Blank 

13-36 Ua6 alpha Curve L Data- bbbbbbbbbbbbbbb 

where b denotes^ a blank 



Columns 

Format 

Nomenclature 

Description 

Card 2 (Curve Header 

Card) 


4-5 

15 

KCRV 

Kind of curve code. Two curves 
may he the same number if kind 
of cvirve is different. 


0 


K curve for pressure drop, 
dimensionless = f(Re x 10“^) 


1 


Density of fluid, or liquid, 

lb /ft3 = f(°F). 


2 • 


(lb /ft-sec)(x 10^) = (°F) 
in 

.. ' '■ ' ■■ 

" 3 • 


Fricti-on factor for fluid', 
f X 103 = f(Re X 10"^) . 
(Used when Re ,2000). 


4 ' ; 


Conductivity, BTU/hr-ft-°F = f(°F). 


5.-' 


Specific heat, .BTU/lb -°F = f(°F) 

m 


6 


Absorptivity, dimensionless = f(°F) 


T 


Emis'sivity, dimensionless = f(°F) 


9 


incident heatj BTU/hr-ft'" = r(hours) 


10 


Prescribed temperature, °F = f (hours 


11 


Total, flow rate, lb /hr = f( hours) 

m 


12 


Fluid inlet temperature, °F = f(hour 


13 


■This card signals the END OF CURVE 
DATA. 

6-10 

15, 

NC 

Curve niunber. 

11-15 

15- 

NP 

Number of points on curve. 

16-72 



May be used for curve title. 

Cards 2 

through 2 NP/7 (Curve Data Cards) 


1-10 

FP 

F10.5 

XI 

Independent variable. 

11-20 


X2 
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Columns Format. Nomenclature Description 

21-30 . X3 

Etc . 

FP Yl Dependent variable. 

F10.5 

Y2. 

, Y3 ' , 

Etc . 

Start Yl in the first field after 
Do not vrite beyond Column 70. 

If the niimber of points given is 1, the value in Columns 11-20 will be 
used for dependent variable. 

U.U.2 Automatic Nodal Subdivision (See Section 4. ^»l) Cards 

The Automatic Nodal Subdivision (ANS) cards are added if NORAD, 
the number of radiator panels to be subdivided, (Parameter Card 2, Columns 
65-66) is not zero. One set of ANS cards are added for each' radiator panel 
to be subdivided. That is, if NORAD = 1, one set of ANS cards are needed; 
if NORAD = 2, two sets of ANS cards are. needed. There are two sets of ANS 
cards, one for rectangular panels and one for circular panels. Each panel 
to be subdivided must be represented by the proper set of cards; ANS-R 
cards for rectangular panels and ANS-C for circular panels. The location 
of the ANS cards relative to the other data cards depends on the value of 
INDATA (Parameter Card. 2, Columns 67-68) which indicates where the data is 
to be found and whether it will be edited by the EDIT subroutine. The ANS 
cards for' side 1 are placed first followed by the ANS cards for side 2. 

The following tabulation describes the location of the ANS cards relative 
to the other data for the various values of INDATA and output tape units 
required for each value' of INDATA. 

LOCATION OF ANS CARDS 
INDATA . 1 THRU 6 


0 After Parameter Card U 

1 After 13 cards 

2 . After last EDIT cards 

(99999 Card) 

3 thru 5 After Parameter Card 2 


OUTPUT 
TAPE UNITS 
REQUIRED 

B 

B.& C 
• B & C 

B 



A description of the MS cards for rectangrolar and circular panels are 
given below. 

Rectangular Panels 

For all rectangular panels (NTYPE(N) = 1, Parameter Card U) to be sub- 
divided by the Automatic Nodal Subdivision Routine, six cards are necessary 


as described below. 



Columns 

Format 

Nomenclature 

Description 

Card MS- 

■ik 



1-10 

E10.3 

X 

Panel length, (x) in. 

11-20 ,N 

. E10.3 

Y 

Panel height, (y) in. 

21-30 

E10.3 

T. 

Panel thickness, in. 

31-33 

13 

NX 

Number of nodes in the x-direction 

34-35 

12 

NTUBE 

Number of tubes in the y-direction 

36-37 

12 

NY , 

Number of . structure lumps between 
tubes. 

38 

11 

CODE 

= 0, parallel flow 
= 1, serpentine flow 

39-52 



Blank - Starting- lump , type and tube 
numbers furnished by TPR. 

53-54 

12 

NZX 

Number of incident heat zones in 
the x-direction. 

55-56 

12 

NZY 

Nimiber of incident heat zones in 
the y-direction. 

57-63 

E7.3 

’HCSk* 

2 

Tube cross-sectional area, in . . 

64-70 

E7.3 

TUBWID* 

Width of tube l\mips in the y- 
direction, in. 

71-76 

a6 

TINIT 

Initial temperature of all lumps , °F, 

77 

11 

IN 2D 

= 1 , 2-D panel . 
= 0, 1-D panel. 


^ Input of both these variables allows the user to distinguish between 
physical tube lump width (TUBWID) and the corresponding dross-sectional 
area. 


10l+ 


Columns 

Format 

Nomenclkture 

•Description 

78 

11 ... ■ 

IN2DC . 

= 1, Connect to following panel 
=2, Connect to first of this series 




of. 1* panels . 

T9 

11 ■ 

IN2DD . 

.“1, Disconnect this panel on the 




diagonal. 

80 



Blank 

Card MS- 

'2R (.Information 

for Fluid Type 

Card) 

1-25 



Blank - Fluid property curve numbers 
supplied by TPR. 

26-30 

5A6 

NKPDC ■ 

Dynamic headloss factor (K) number. 

31-50 

Blank 



51-60 


CSA .. 

2 

Fluid cross-sectional area, in . 

61-67 

3A6 , Ah 

WP ' 

Wetted perimeter, in. 

68-72 


FRE 

Friction factor adjusting, co- 
. efficient. 

Card MS- 

•3R (Information for .Tube Tyj>e 

Cards ) 

1-5 


DENM 

Density, Ib^/ft. 

6-10 


NCONC 

. Conductivity curve number. 

11-15 

ltA6, A1 

. NSHC 

Specific heat curve nxunber. 

16-20 


NABC 

Absorptivity curve number. 

21-25 


NEMC, 

. Emissivity curve number. 

39-itO 

12 

LCC 

Longitudinal conduction code 
= 6, Ivimps in this panel conduct 




longitudinally 




= 1, lumps in: this, panel, do not 




conduct longitudinally 

Card ' MS- 

•4r (More information for, tube 

type cards) 

1-10 

3A6, A2 

MT 

Area for heat transfer to 
enclosed fluid lump, in^. 

li-20 . 


. AE . ■ 

Area of surface for external : 
. radiator, in^. . 

21-30 

E10.3 

YCTS ■ 

Conduction distance for tube 


to structure, in. 



ColiAmns . Format 


Nomenclature 


Description . 


Card ANS-3R (information for Structure Type Cards). 

I- 5 DEN Density, 

6-10 NCONC Conductivity curve number 

II- 15 Ha 6, A1 . NSHC Specific heat ’ curve number 

l6-20 NABC Absorptivity curve number 

21-25 NEMC Emissivity curve number 

Card ANS-6R 

1-5 .15 Incident heat curve number for 

' - . * - • ' - "- -first zone. - ” 

6-10 15 . Incident heat curve number for 

second zone. 

Etc .to 

1*6- 50- 1 5 Incident heat curve number for 

tenth zone. 

Repeat Card ANS-6R as many times as necessary to supply an incident 
heat curve number . for each zone'. 

Repeat Card ANS-IR through Card ANS -6R for each rectangular radiator 
panel to be subdivided. 

Circular Panels 

For each circular panel (NTYPE(N) - 2, Parameter Card i*). to be subdivided 
by the Automatic Nodal Subdivision Routine,' six cards are necessary as 
described below. 

Card ANS-IC 


1-10 

FIO . 0 

Radius 

Panel radius, in. 

11-15 

F5,.0 

DELTA 

Panel thickness, in. 

16-20 

F5.0 

DIAMI 

Inside tube dimeter, in. 

21-25 

F5.0 

LMAX 

Maximum length of tube node 
desired, in. (may be left blank) 

26-30 

F5.0 

LMIN 

Miiiimum length of tube node allowed 
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Coliimns 

Format 

Nomenclature 

Description 

36-37 

12 , 

NT 

Number of tubes on panel. 

38-52 



. Blank - Starting lump , type and 
tube numbers furnished by TPR. 

,53-54 

12 

NS 

Number of structure lumps between 
tubes. 

55-56 

12 . 

NQL 

Incident heat curve number for last 
aone. 

57-63 

F7.0 

TCSA 

Tube cross-sectional area, in^. 

64-70 

>7.0 

, TW . 

Width of tube .lumps in radial 
direction, in. 

71-79 

F9.0 

TINIT-. 

Initial temperature of all lumps , °F 

80 



Blank 

Card MS- 

2C (Needed only 

if LMAX < 0.0) 


1-5 


NTBS(l) 

Number of .tube lumps for first tube. 

6-10 


NTBS(2). 

Number of tube; lumps for secohd 
tube. 

11-15 

' 15 . 

, NTBS(3) : 

Number bf tube liamps for third 
tube. 

Etc .to 




76-80 

15 

NTBSC16) . 

Number of tube lumps for sixteenth 
tube. 

Repeat card MS-2C as many times as necessary to supply number . of tubes 
lumps for each tube. 

Card MS- 

3C (Information for Flbid Type Card) 

1-25 



Blank - Fluid property curve numbers 
supplied by TPR. 

26-30 

I5'. 

’ NKPDC :: 

Dynamic head loss (K) factor n\imber. 

31-50 



Blank 

51-60 

FIO.O 

CSA 

\ . '■ .2 . ■ ' 

Fluid cross-sectional area, in . 


lOT ■- •; . 



Columns 

Format • 

Nomenclature 

6i-6t . 

FT.O 

WP 

60-7;.' 

I'"] .0 

FHE 


Card MS-Uc (information for Tube Type 


1-5 

F5.0 

DENMT 

6-10 

15 

NCONCT 

11-15 

15 

NSHCT 

16-20 


NABCT 

21-25 

-15 - -- 

= --NEMCT 

39-^0 

12 

LCC 


Card ANS- 

5C (Information 

for Structure 

1-5 

F5.0 

; DENS 

6-10 

15 

NCONCS 

11-15 

15 

NSHCS 

16-20 

15 

NABCS 

21-25 

15 

MENCS 

Card ANS- 



1-10 

FIO.O 

YICEPl 

11-20 

FIO.O 

XICEPl 

21-30 

, FIO.O 

BETAl 

31-35 

■ 15. .. 

NQCURV 


Description 

Wetted perimeter.., 

Frio lion factor adjusting co- 
efficient. 

Cards ) 

Density, lb /ft^ 
m 

Conductivity curve number. 

Specific heat curve number. 

Absorptivity curve number. 

- - Emissivi-ty curve- number. 

Longitudinal conduction code 
= 0, lumps, in this panel conduct 
longitudinally 

= 1, lumps ■ in this panel do not 
conduct longitudinally 

Type . Cards ) 

Density, lb /ft^ 

Conductivity curve number. 

Specific heat curve number. 

Absorptivity curve number. 

Emissivity ctirve number. 

Y-axis intercept for incident 
heat zone boundary line number 1. 

X-axis intercept for incident 
heat zone boundary line number 1. 

Angle between incident heat zone 
boundary line and (+) x axis of 
pane], degrees, (Needed only if 
YICEPl = XICEPl = 0.0) 

Incident heat zone curve number. 
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Colijmns 

Format 

Nomenclature 

Description 

36-40 

15 

NLINE 

Number of boundary lines defining 
, incident heat zone (1 or 2). 

4i-50. 

■FIO.O 

YICEP2* 

Y-axis intercept for incident heat 
zone boundary line number 2. 

51-60 

FIO.O 

XICEP2* 

X-axis intercept, for incident heat 
zone boundary line number 2. 

.61-70 

FIO.O 

BETA2* 

Boundary line angle (degrees). 

75 

11 

NCODE 

= 0, Indicates more incident heat 
zone i nfomat ion follows . 

= 1, Indicates last incident heat 
zone data card, 

= 2, Indicates only one incident 
heat zone for this panel. 


Repeat Cards ANS-6C as many times as necessary to supply incident heat zone 
data for each zone. 

Repeat - Cards ANS-IC through MS-^C for each circular radiator panel to be 
subdivided. 


* These parameters are needed only when NLINE = 2. 



Two-Dimensional Rectangular Panels 

For each 2-d panel (NTYPE (N) = 3» Parameter Card, U) to be 
subdivided by the Automatic Nodal Subdivision Routine, five cards are necessary 
as described below. 



Card ANS-ID 



Columns 

Format 

Nomenclature 

Description 

1-10 

FlOlO 

PANWID 

Panel width, in. 

11-20 

FlOlO 

TUBWID 

Tube width, in. 

21-25 

F510 

DELTA 

Fin thickness, in. 

26-30 

F510 

TCSA 

2 

Tube cross-sectional area, in. 

31-35 

F510 

DIAM 

Inside tube diameter, in. 

36 -Uo 

F510 

TSPl 

Tube spacing for first tube (leave 
blank for equal tube spacing) 

41-45 

F510 

XMAXl 

Maximum length of tube nodes for 
first tube (may be left blank) 

46-50 

F510 

XMAX2 

Maximum length of tube nodes in 
remaining tubes (may be left blank) 

51-60 

FlOlO 

TINIT 

Initial temperature of all lumps 

61-65 

15 

NT 

Number of tubes on panel 

66-70 


NTia 

Number of triangles on panel 
(Must equal 4) 

71-75 

15 . 

NDIUA 

Number of allowable subsequent 
tube divisions after first tube 
(overrides XMAX2 option, may be. 
left blank) 

76-80 

15 

NQZONE 

Incident heat zone number 


Card ANS-2D (Use only when XMAX2 =0.0) 

1-5 

15 

NTLl 

Number of tube lumps in first tube. 
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Columns 

Format 

Nomenclature 

6-10 


NTL2 

76-80 

15 

NTLI6 


Continue 

on another 

card if NT 16. 

Card ANS-3D 


1-5 

15 

HKPDC 

6-10 

11-20 

FlOlO 

FCSA 

21-30 

FlOlO 

WP 

31-^40 

FlOlO 

FRE 

Card ANS-4D (Tube type data) 

1-10 

FlplO 

DENT 

llrl5 

15 

KT 

16-20 

15 

CPT 

21-25 

15 

ALPHAT 

26-30 

15 

EMIST 

31-35 

15 

LCC 


Description 

Number of tube lumps in second tube 

Number of tube lumps in sixteenth 
tube. 


Curve number for N^’PDC curve 
(May be left blank) 

Blank 

. ' ■ ‘ ' 2 

Fluid cross-sectional area, in. . 

Wetted perimeter, in. 

Friction factor adjusting 
coefficient (nor^lly = l.O) 

Density of tube lumps, Ib^/ft 

Conductivity curve number 

Specific heat curve number 

Absorptivity curve number 

Emlssivity curve number 

longitudinal conduction code 
= 0 considered 
= 1 not considered . 
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Card ANS-3D (Structure type data) 


Columns 

Format 

Nomenclature 

Description 

1-10 

FlOlO 

DENS 

Density of structure lumps, 

lb /ft^ 
m' 

U-15 

15 

KS 

Conductivity curve number 

16-20 

15 

CPS 

Specific heat curve number 

21-25 

15 

ALPHAS 

Absorptivity curve number 

26-30 

15 

EMISS 

Emissivity curve number 


U.U.3 Plot Program 

1. Control Cards 

. T/8 ASG A = ,PCF (Program tape) 

7/8 ASG I = PLOT (History tape). 

.. 7/8 XQT CUR 

TRW A 
IN A 
TRI A ■ 

7/8 XQT PROG 

Parameter cards and data cards, if any 
7/8 XQT PLOT A . 

Plot data card's 
7/8 EOF 

2. Data Cards 

Columns Format Name Description 


a. Case 

Title Card 



1-72 

12A6 

■:TITLEA 

Title to be printed at .the top of each 




grid. 

b . Time 

Card 



1-10 

■ FIO.O 

"TA 

First value of time to be plotted (hrs) 

11.- 20 

FIO.O 

TZ . 

Last value of time to be plotted (hrs) 

21-30 ' 

FIO.O 

. . TPG 

Time range for each grid (hrs) 

c . Items 

1 Cards (Any quantity. 

followed by a blank card) 

1-5 

15 

. ITEM 

The item niimber to be plotted. Use a 


negative value if this item is to start 
a new grid. A maximum of four curves 
may be plotted on one grid. Insert a 
blank card when the number of items exceeds 
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Columns 

Format 

Name 

Description 




35,000 divided by tlie number of points 
between TA and TZ. More, item cards may 
then follow this blank card if more curves 
are to be plotted. 


A2 . 

, ITYPE 

A two character item type code. (See 3 
for definitions ) 

11-58 

8a6 

TITLES 

Item description to be printed at the top 
of each grid, along with the plotting 
symbol which is generated and used by the 
program. 


The next two values. are optional on. the cards whose item numbers are negative 
and are ignored on all other item cards. 


6i-to 

FIO.O 

■ YLO 

The minimum value of the y-axis . 

71-80 

FIO . 0 

YHI 

The maximum value on the y-axis . 


3. Detail Type and 

item Descriptions 

Type 

Item 


Description 

MP 

1 


Pressure drop in system 

PR 

1-N 


Pressure drop of 1 through N tubes 

BT 

1 


Heat rejection of system 1 

FR 

1-N 


Flowrate of 1 through N tubes 

FT 

1-N . 


Fluid temperature of 1 through N lumps 

TT 

1-N 


Tube temperature of 1 through N lumps 

ST 

1-N 


Structure temperature 1 - through N lumps 


llU 



5.0 


COMPUTER INCIDENT HEAT ROUTINE 


The Computer Incident Heat Data Routine (CIHR) was developed 
to reduce the large amount of user time required to generate incident 
heats for use in LTV thermal analyzers for orbital missions . Prior 
to the availability of this routine this task consisted of (l) running the 
Midwest Research Institute (MRl) Routine to obtain the incident heats, (2) 
selecting the number. of points required to describe the curves, and (3) 
filling out data sheets for key punching the curves in a format compatible 
with LTV routines. Using the CIHR, only the first step is required. The 
remaining steps required to obtain the incident curves on cards are per- 
formed by the computer including (l) siirveying the points to determine the 
minimum number needed to define the c’urves and (2) punching the cards in a 
format compatible to LTV routines . 

The CIHR was created by modifying the Midwest Research 
Institute’s Heat Rate Routine described in" Reference U to : (iL) survey 

the data points of each heat rate curve generated by it and select only 
the points necessary to describe the curve for linear interpolation, (2) 
punch an output curve generated from the selected points on cards with 
the appropriate curve number in LTV curve format, and (3) plot the curve 
created by the selected points with original curve on the same grid and 
the output curve on another grid. The method .used to. select points along 
with the additions to the data preparation described in Reference 5 are 
given below. 

Method 


. Points to describe a heat rate curve are selected from the 
points of the curve itself such that the heat rate difference between the 
original curve and curve created by the selected points is never greater 
than a tolerance. In Figure 5-1 » points A, B, C, and D are the data 
points of the original curve. Points A, C, and D are selected to describe 
the curve. Point B is eliminated because Dl. is less than the tolerance 
and C is selected because D2 is greater than the tolerance. 

Data Preparation 

The data preparation for the computer incident heat routine 
is given in Reference 5 except for the following new values on the low 
altitude, planet temperature card and element card. 

A list is given on the following pages summarizing the; 
various data that might be punched as program input data. For each type 
of card a description is given of all applicable fields. Page niombers 
referencing additional information in this manual are also given. If the 
entry under "Applicability" is blank, the field in question may pertain 
to all kinds of cases. If there is an entry under "Applicability", however, 
the field is pertinent only for the type of case designated. For all other 
types of cases the field is not applicable, i.e., it can be left blank. 


115 . 


Time 


FIGURE 5-i ILLUSTRATION'GF METHbD USED TO SEiIeCT POINTS 





Note that the entire 05 card is not applicable to trajectory, tape missions 
and can be omitted in such cases. Likewise the 08 card is not applicable 
to orbital missions. 

The formats listed are those used by the program; however, for 
all fields specified as E8.1 it is advisable to punch the number with a 
decimal point but with no exponent unless the magnitude of the n\amber is so 
large or so small as to require exponents. For example, one (l) in this 
format could be punched as "O.IE 01’*, but fewer errors would be made if it 
were simply punched as "liO" anywhere; in the eight-colxann field. Numbers 
in. I format must be right Justified in their . appropriate, fields . 



Columns 

1-2 

3-5 

9-12 

13-20 

21-28 

37-44 

78 

79: 

8d 


CASE NUMBER AND CONTROL INFORMATION (Ol CARD) 


Ref. 5 


Format 

Applicability 

Description 

Page 

12 


Card code (=01). 


13 


Case number. 


I4 


Number of elements to be analy- 
zed (14n4l000). 


E8.1 

Orbital 

0Q, the initial value for true 
anomaly. 

49 

E8.1 

• Orbital 

0Q, the increment of true 
anomaly at which heat 
fluxes are calculated. 


E8.1 

TraJ . tape 

. First time (min. ) in which pro- 
. gram parameters are to be in- 
stantaneously changed. . 


11 


Configuration number (1-6) if 
Apollo shadow factors stored 
on tape are required for thi? 
case. (Shadow tape must be 
mounted when a valid' configura- 
tion number is used.) Blank 
or 0 is used when no shadow 
factors are used. 


11 


For Oor blank, only absorbed 
heated are printed out. 




For 1, incident and absorbed 
heats are printed out. 


. 11 ' 


For 0 or blank, fluxes in 
BTU/hr-ft^ are printed out. 



For 1, products of flvixes and 
element areas in BTU/hr are 
printed out . . 





. ORIENTATION AND ROTATION (02 CARD) 


Coliimnri 

l'’orma.t 

Applicability 

Description 

Ref. 5 
Page 

1-2 

12 


Card code (=02) . 


3-5 

13 

Rotating 

fl (°), one of angles describ- 
ing axis of vehicle rotation. 

53 

6-8 

13 

Rotating 

A (®), one of angles describ- 
ing 6Lxis of vehicle rotation. 

53 

9-12 

lU 


For 1, vehicle is planet orient- 
ed.. 


. . • . — 

- ; 

... ■. .... 

For 0, vehicle spins rapidly 
about random-axes ■ - 





For -1, vehicle is sun-oriented. 





For -2, vehicle is star-orient- 
ed . 


13-20 

E8.1 


(°), vehicle roll angle. 

52 

21-28 

E8.1 


p(°), vehicle pitch angle. 

52 

29-36 

. E8.1 


'|^(°), vehicle. yaw angle. 

52 

3T-4U 

E8.1 


17 (rph), rate . of vehicle ro- 
tation. Blank or 0 is used 
when vehicle is not rotating. 


^5-52 

E8.1 

Star- 

Oriented 

X-component of a vector pointing 
to orientation star. 


53-60 

E8.1 

Star- 

Oriented 

Y-component of a vector pointing 
to orientation star. 


61-68 

E8.1 

Star- 
Oriented . 

Z-component of a vector pointing 
to orientation star. 



PLANET (03 CARD) 


Columns 

Format 

Applicability 

Description 

1-2 

12 


Card code (=03). 

3-5 

13 


For 1 to 9> number is taken as 
planet code. 

8. 

11 


0 for constant surface temp. 

1 for variable temp. 

For 0 or blank, planet data 
are read from columns 13-52. 

13-20 

E8.1 

Unknown 

planet. 

Distance from. planet to sun in 
nautical miles. 

21-28 

E8.1 

Unknown 

planet 

Planet radius in nautical miles. 

29-36 

E8.1 

Unknown 
planet . 

Planet albedo. 

37-^^^ 

E8 .1 

Unknown 
planet . 

3 2 

GM (ft /sec ), gravitational 
constant times planet mass. 

^ 5-52 

E8.1 

Unknown 

planet 

Adjusted cold side temperature 
(°R). 
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Rel'. 5 
Page 


60 


MISSION (OU CARD) 


Ref. 5 


Columns 

Format 

Applicability 

Description 

Page 

1-2 

12 


Card code (=0U). 


3-5 

13 


For 0 or blank, trajectory tape 
mission is to be run. 





For 1, orbital mission is to be 
run . 


6-8 

.13 

Orbital 

Number of degrees of true anomaly 
for which program computes. 

h9 

9 , ■ 

11 - 

Tra=j . . tape - - 

= Planet code of first^ trajectory 
tape planet of reference.^ 


10 

11 

Traj . tape 

Planet code of second trajectory 
tape planet of reference.* 


13-20 

E8.1 

. Orbital. 

Apogee in nautical miles. 


21-28 

E8.1 

Orbital 

Perigee in nautical miles. 


* . Unless 

new routines. are written 

to manipulate the trajectory tape 

(see 


Appendix I of Reference 3-1 ) the. planet codes are 1 and 2, respectively. 


SUN POSITION (03 CARD) 


Coiumna 1'’orihat Applicability Description 

1-U 12 Orbital . Card code (=05). 

3-5 13 Orbital A code to specify sun position 

input option. Value must be 1, 
. 2, 3 or 5- The code instructs 
the program how to interpret 
data in 13-T6. 


Ref.. 5 
Page 


6-8 13 Orbital Fprl, sun-shade points 0. and 

0 . are read from columns o 1-T6. 
Otherwise they are computed. 

13-20 . E8.1 Orbital • Angle i(°) if column 5 contains 51 

1 or 5. 

'Angle5^(°) if column 5 contains 50 

2 ,. 

Angle .(x{°) if coliunn 5 contains 50 

,•■3.' , ■ 

21-28 E8.1 ■ Orbital Angle (u (° ) If .column 5 contains 51 

1. dr 5.- ' 

Angle /3 (°) if; column 5 contains 50 

2 or 3. . ■' 

29-36 E8.1- Orbital . Angle.- n (®) if colimn 5 contains 51 

1 or 5;. 

Angle y (®) if column . 5 cpntains 50 
3. Otherwise blank. 

37-44 E8.1 Orbital Angle RA(®). if column 5 contains 50 

Year if column 5 contains 1. 

Otherwise blank. 

45-!-52 E8.i Orbital ^ Angle DEC(®) if column 5 contains 50 

Month if column 5 contains 1. 

Otherwise blank. 
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Columns 

53-60 

61-68 

69-76 


. SUM POSITIOM (05 CARD) (Concluded) ' 

Ref. 5 

Format Applicability Description Pa^e 

^8.1 Orbital. Day of month if column . 5 contains 

1. , Otherwise blank. 

E8.1 Orbital Hour. of day (O.to 23) if column 

5 contains 5- 

0. (°) if column 8 contains 1. 

O^Serwise blank. 

E8.1 Orbital Minute of hoiir (0.. to 60 .) if : 

' column 5 contains 5 * 


' — 0 r ( ° )' if" column* 8- contains 1 . 
Otherwise blank. 


Columns 

1-2 

3-5 

6-8 

9-12 


13-76 

77-78 

1-2 

3-5 

6-8 

9-12 


REDEFINITION (06 CARD) 

Ref. 5 

P’ormat Applicability Description Page 

(Subcase with 3, 5» or , 6 in coliamn l2) 

12 Caird code (=06). 

13 Index of first table entry on card. 

13 . Index of last table entry on card. 

lU , Code describing type of table being 

defined. 

Code 3 indicates f versus t. 

Code 4 indicates f versus t. 

Code 5 indicates g versus 8 . 

Code 6 indicates b versus T. 

m. , 

8e 8.1 Table entries to be stored as index 

ranges from first to last value. 

12 Number of table being defined 

(from 1 to l6) . 

(Subcase with 7j 8, or 9 in coliamn 12) : 


12 

Card code (=06). 

13 

Material code (from 1 to l6). 

13 

Corresponding table numbers 
. (from 1 to l6) . 

l4 

.Code describing type of table 
assignment. 


Code . 7 indicates assignment of 

f versus t table. 

P 


Code 8 indicates assignment of 

f versus t table, 
s 


Code 9 indicates assignment of 
g versus 8 table. 
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REDEFINITION (06 CARD) (Concluded) 

Ref. 5 

Columns Format . Applicaljility Description Page 

(Uuljoase wiLfi 10 in columns 11 - 10) 


1-2 

12 

Card code (=06). 

9-12 

lU 

Subcode (=10) indicating solar 
constant is redefined. 

13-20 

e8.1 

New value of solar constant. 
(Subcase with 11 in columns 11 - 12) 

1-2 

12 

Card code (=06). 

3-5 

13 

Code of first material referenced 
on this card. 

6-8 

13 

Code of last material referenced 
on this card. . 

9-12 

lit 

Subcode (=11) indicating that this 
card redefines values of a for 
one or more materials. ^ 

13-T6 

8e8.i 

Values of a for materials with 
codes ranging from first to last 


as shovn in columns 3-5 and 

6 - 8 . 
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Coilunns 

1-2 

3^5 

9-12 

13-20 


21-28 

29-36 

37-i+^ 

> 5-52 

53-60 

61-68 

77-78 

79 

80 


LOW ALTITUDE PLANET TEMPERATURE (07 CARD) 

Format Applicability Description 

(or Name) 

12 Card code (=07). 

13 NSELC =0, No output curves requested. 

=1, Output curves requested for 
specified elements. 

iL NTOL . lype of tolerance code for CIHR 

=0, Absolute. 

= 1 , Relative. 

E8.1 ■ Altitude (n.m.) below which planet 

temperat\ires computed by, the pro- 
gram are superseded by temperature 
■ in columns- 21-28 of this card. 

For blank or 0 . the program always 
uses computed temperatures. This 
altitude is understood to be zero 
. unless an 07 card is read to re- 
define it. 


E8.1 


Planet temperature (°R) used by the 
program when the altitude is, below 
the value in columns 13-20 of this 
card. 

E8.1 

TOL 

Tolerance for CIHR. (BTY/hr if NT0L= 
0 and percent if NTOL =1). 

E8.1 

START 

Start time for output curves (hr) 

E,8.1 

STOP 

Stop time for output curves (hr) 

E8.1 

. TiMADD 

Time adjustment to the; original 
curve (hr) 

E8.1 

QADD 

Heat rate adjustment to the ori- 
ginal c\irve (BTU/hr). 

12 

NORBIT 

Number of cycles of selected data 
. points to be in output curves . 

11 . 

NPUNCH 

= 0 , No, punched output 
= 1 , Punch output curve on cards 

ii 

NPLOT 

= 0 , No plot output 
= 1, Plot output 


Ref. 5 
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If no 07 card is read, no output curves are K^nerated. If NMlilliC oqun.'li' 
zero, the other variables described above are not used, if START is 
less than the time of the first data point, the. output curves have zero 
heat rate from START to the first data point. The output curves end 
at STOP or after NORBIT cycles of the selected data points whichever 
is greater.. If STOP is greater than the time of the last data point, the 
selected data is cycled until STOP is reached. 


TRAJECTORY CONTINUATION (08 CARD) 


Columns Format . Applicability . Description 


Ref. 5 
Page 


1-2 12 Traj , tape Card code (=08). 


13-20 E8.1 Tra j . tape Cutoff time (min.) for next seg- 

ment of trajectory. When a time 
^ ‘ ‘ - exceeding this value is read -from 

the trajectory tape, the case . 
will be interrupted to read para- 
meters being instantaneously 
changed. 
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. ELEMENT (09 CARD) 

Ref. 5 

Columns Kor*mat Applicability Description 

1-2 12 Card code (=09). 

3-5 .13 Sequence number of element for 

this problem. The momber must 
be no greater than the number of 
elements to be analyzed. Elements 
. to be lumped together require 

adjacent sequence numbers. 

6-8 . 13 If any of the internally stored 

element data are needed’, the node 
number for the element must be 
punched. If this field does not 
contain a valid node nxjmber area, 

A and ft are read from this 
''' cardj and sRadow factors cannot 

be referenced. 

9-12 ih' This field contains a nonzero 

numerical label which will identify 
the element in the output listing. 

If this field is blank or 0 the 
program uses node- number in 
columns 6-8 as label. 

13-^20 E8.1 Angle A (®) for the element in ^9 

case values from internal tables 
are not used. 

21-28 E8.1 Angle ft for the element in ^9 

case values from internal tables 
are not used. ' ' 

29-36 E8.1 Element area (ft ) in case values 

from internal tables are not used. 

37-^^ E8.1 For all elements which are to be 

lumped with succeeding elements, 

. i.e. , with elements having succeed- 
ing sequence numbers, this field 
must contain the same label punched 
in columns 9-12 unless they are 
blank, in which case it must con- 
tain the node number from columns 
■ 6-8. Note that the format here 

requires a decimal ppint. 
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. Columns 



53-60 


61-68 


69-76 

TT-T8 

T9 


80 


. ELEMENT (09 CARD) (Concluded) 

Format /Vpp I inability Description 

E8.1 For 1.0 the value of A in 13-^0 

supersedes the value stored in in- 
ternal tables. . For blank or 0. 
the value, from the internal table is 
used if a valid node number is 
given. 

E8.1 For 1.0 the value of Q in 21-28 

supersedes the value stored in in- 
ternal tables. For blank or 0. the 
value from the internal table is 
used if a valid node number is 

E8.1 For 1.0 the value of area in 29-36 

supersedes the value stored inter- 
nally. For. blank or 0. the value 
from internal tables is used if a 
valid node number is -given. 

E8.1 NCURVE 0, Curve n\amber of output curve for 

this element. 

0., No output curve for this element. 
Note that the format here requires 
a decimal . point , 

12 Coating material number. 

II For 1, shadow factors . for this node 

are loaded from the shadow tape and 
shadow logic is used to find heat 
. fluxes. For 0, shadow logic is 
. suppressed and factors are not 
loaded. 

II For 1, heats will be plotted for 

this element. For 0, plots will be 
suppressed. 


Hof. ‘i. 
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. COMMENT (10 CARP) 

Coliirnhs Format Applicabiiity . Description 

1-2 : 12 . . Card code (=10), 

3-80 13A6 Alphameric . comment . 


Ref. 5 
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APPENDIX A 


DETERMINATION OF RADIAL FIN EFFECTIVENESS 
IN A NON- ZERO SINK TEMPERATURE ENVIRONMENT 



1.0 


SUMMARY AND INTRODUCTION 

This Appendix presents the method and results of an analysis 
which was conducted to determine the fin effectiveness for circular radiating 
fins in a non-zero sink temperature environment. The second order nonlinear 
differential equation which describes the heat flow within the fin along with 
the associated boundary conditions were derived from the physical laws of 
heat transfer.. The solution of the governing equation had been determined 
prior to this study for a limited range of design parameters, (Ref. lA ), but 
had not previously been determined for non-zero sink temperatures. In this 
study a numerical teclinique was devised whereby values of radia.l fin ef- 
fectiveness could be determined from the governing equation for sink tempera- 
tures greater than or equal to zero. A computer program. Radial Fin Effec- 
tiveness Routine (RFER), was written incorporating this numerical technique 
and the solution was obtained for a wide range of design parameters with the 
use. of a NASA-MSC Univac 110 8 computer. A detailed discussion of the method 
and results is presented below. 

2.0 DISCUSSION 


The governing differential equation for a circular radiating fin 
can be derived.as follows. Consider a heat balance on the surfaces of the 
differential element of a radial fin shown in Figure A -1, 


On surface 1: 

5t 

Q = -kA ^ 

Q = K 5 rdQ 
r 0 r 


A = 5rd0 


On surface 2: 


Q . . = Q t -r — (Q ) dr 

r+dr r • dr r 


where: 


On surface 3: 

«rad = “"'r ' 

Q , = Ns f ardO dr (T^-T„^) 
raci o 


Q = heat flux 

Ns = mamber of sides radiating 
^ = thermal emissivity 
5 = fin thickness 
K = thermal conductivity 
o = Stef fan-Boltzman constants 


A = rdOdr-Ns 
r 


A-2 




Assuming one dimens'ioriai, radial heat flow, a Heat balance gives: 


“:l„ = “cut 


% ' “r+dr * ®rad 


Q ^ - Q + Q , = 0 . 

r+dr r rad 

dr +fard9dr (T— T^)Ns = 0 
dr . S 

[-KS rdQ -^1 dr + f ard0dr(T^- Tg^)Ns = 0 

-K 5 2 rr-^ [^ 4^] dr +. 2 rreardr (T^ - T„^) Ns =. 0 
a r a r S 

-K S [ If , * . +«7r(T*' - Tg‘‘)Ns = 0 


^ . 1 _M . T**) =0 

ar^ r or oK S 


Eq. 1 


If we let: 

1) 0 = T/Ti 

2) R = 

ro-ri 

3) 0s = Ts/Ti 

1.) <r = Tl3(r^.rj)2 

K S 


Ti = Temperature at r =' ri 


5) y = - 


r _ r 
b i 


6) 0 = 


O 


r . 

1 


r = ri + R(r - r. ) 

o 1 

h h h 

T = 0 T. 

1 • • 

. h h k 
^S.= V^i , 


(T^ T = T^^(0^ - 0^^). 
S .1 s 


dr 


'T^d(T/T^) 
(rp-ri) d,r-r. , 

. r -r . 
o 1 


T. . 

<v^i> 


dR 


A-U 



ill 

dr 


d^T _ 


' 1 ’,. 

(r -‘ r. ) 

U I 

T. 


dr 


d9 

dU 


r^l = i 

•- ATD J 


d(dO) 


(r - r. ) dr ^dR-^ (r -r. ) (r -r. )d(r-r^ ) [dR] 
o 1 oi oi i- 


r -r. 

o 1 


o 

d' T 

■■■ "O’ 

dr"' 


T. 

1 




d^e 


dR 


equation 1 gives : 


Substituting the above values of , r, (T^-T and d^T into 

dr dr2 


- T - ■■ • ™ - - - - -■ - - T* ~ - - - r— ' I*'-/ i* ■ . 

i ^ p- ■ 1 . i ^ _ Ns e crT. (9 -9s ) 

/ T2~ dr^ r. + R(r -r. ) / x dR jTT 

(r -r. ) 1 oi{r-r.) K5 

O 1 O 1 


= 0 


d^9 

— — n’ + 

dR^- 


(R .+ : r. /r -r. ) 
10 1 


^ NseaT^^(r -r. )^(9^-9s^) 

■' ■ ' .Ks° " 


^ -H -1. - C (Q^-Qs^) - 0 

d^ (R + y) dR ■ ® ° 


Eq. 2 


Equation 2 is the governing differential equation, for the tempera-^ 
ture of a one-dimensional radial fin in dimensionless form. 


Calculation Of. Q from fin: ■ 

Q= rjfa^sn (r^^-r^^) (T^^-Tg^) 
Calculation of Q conducted from tube: 


Eq. 3 


Q = -KA (^) 
dr 


r=r. 

.1 


= -K.a 2 . r.,(f) 


r=r. 


•1 


Eq. h 


Setting T. = 

1 13 


9 = 




R= 1:21. 

r -r. 
o 1 


dr 


* dR “ / \ 

<Vh> 


dr 


r . 

1 


• 1 


A-5 



dO dT 


dl<' dr T 


h 


B 


(r -r. ) ''dR^ 

r=r. o 1 R=0 

1 ■ - 


(iM) 


Setting equation 3 = equation h 

r, 2 2. k ^ h. -K 5 2 77 r.T^ ' .dQ. 

7}€0]is n (r -r ) (T -T_ ) = i B (— ) 

. O 1 o o / \ Un 

(r -r. ) R=0 

o 1 


/dQy _ - 17 f <7 Ns 77- (r ^-r ^)(r -r )(T 

V jp J “ ■ o i o i B S 

. R=0 K S 27rr,T^ 


i B 


(^) 


R=0 


neaNsir -r. )^T^^(l-9„^) (r +r. ) 
. ■ ■ 01 B S o 1 

2KSr. 


(^) 

^dR'^ 


77 C d/2 + V )(l-0g^) 


R=0 


■ Rq. S 


l'’roin i'jq . : 

2 

/d 0. h 

= <d-eg) 

R=0 ® 


(^) 


(o+y) 


( 4 ) 


C(i - o. - 


R=0 


(^) 


R=0 


Eq . 6 


P 2 

Now with the values of d0/dR and d 0/dR at the ‘ fin Base, (R=0) 
known the differential equation can he integrated in step-wise manner. The 
only parameter which is not known is ij or fin effectiveness,. The procedure for 
obtaining the correct value of 7j is outlined in following section. 
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3.0 


equation 


PROCEDURE FOR CALCULATING RADIAL FIN EFFECTIVENESS 


Governing differential equation in dimensionless form: 


n 

(i^O 


I- i!8 / e S - n 

dR^ (R+y) <i.R C -63 ) - 0 


Eq. T 


Boundary Conditions: 

1 ) 9 ( 0 ) = 1 

2) O' (1) = 0 

3) 0 (1) > 0 


Radial fin effectiveness is obtained by a step-wise' integration of 


Calculate for each R. ; R. = 0, 1, 

1 1 ’ 


®i^l = ®i ^ (f)i 


Eq. 8 


(^) = (^) 
i+1 ^ i 


,d^9. 


car) 


Eq. 9 


\ = d: fo^-o 1 

\ ^ . s' (R.^n + y ) 


i+l 


i+1 


(^) 

^dR\ 

1 


Eq. 10 


Calculate .constants and 6 0 " 

. 0 o 


1 ) 

2 ) 


y 


3) (— ) 


R=0 




R=0 


5) Make an initial guess of ri and go through equations 8-10. 

6) Check dO/dR against zero at the end of the fin (R=l). 

7 ) If not within small tolerance decrease rj and repeat process 
until dO/dR approaches zero. 
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One must be very careful in arriving at the initial guess of rj . 

If initial 77 is chosen below the actual value, then of course the process 
will not converge on the correct answer. 

Tliri'i' arc* l‘t>ur VM. 1 ol’ rj Ui.mI. will ylcM m. .*il. l\ 1. 

Tlio correol. va.l.ue of f In erfectiveness will satisfy boundary condition 3 
while the other 77*5 yield negative values of 9 at the end of the fin. 

A computer routine was written capable of computing values of radial 
fin effectiveness for nearly the entire range of design parameters. Tlie numerical 
•technique described above was programmed into the Radial Fin Effectiveness 
Routine (RFER) and used to calculate values of 17 over the range shown below: 

1) = 0., .5, . 6 , . 7 , . 8 , .9 

2 ) ^ = 0 . 1 - 200.0 

. 3 ) = 1.5, 2.0, 3.0, 6.0, 10.0, 15.0, 50.0 

Curves of the results are shown in Figiires A2-A13- 

The curves of radial fin effectiveness were curve fit for addition 
into SSDR. The curve fit was based on a modification of the value of 77 of 
a one-dimensional fin at the stme value of ^ and 0^. Tills method fu-ovider* .*i 
simple curve fit since SSDR already had 1-D fin effectiveness c;ipab;i In ty . 

0 < ,^ < 1.0 X 

1.0 < C < 4.0 X 

k.o < C < 15.0 ■ X 

15.0 <^< 200.0 / X 

^ - fill effectiveness for a 1-D fin 

X 

V = ^ ~ Curve fit equation programmed into SSDR. 

ri 

U.O . REFERENCES . ' 

1-A Ch^bers, R. L. , and Somers, E. V. , "Radiation Fin Efficiency For 

One-Dimensional Heat Flow in a Circular Fin", Journal of Heat Transfer , 
November 1959. 


= -. 0 T 2 i+ In { ( ) - . 2)^0 
= -.075 In ( C ) - -248 
= -.OU786 In ( (^ ) - .215 
= -.0539 in ( C ) - .33U 
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APP^IX B _ 

DETERMINATION OF FIN EFFECTIVENESS FOR 
A HEAT ABSORBING FIN 


1.0 


SUMMARY AND INTRODUCTION 


This appendix presents the method and results of an analysis 
which was conducted to determine values of fin effectiveness for an 
absorbing fin. The differential equation describing the heat flow in a 
one-dimensional radiating fin was integrated for a reasonable range of 
boundary conditions to obtain effectiveness for values of sink temperature 
greater than the fin base temperature. A curve fit of the information was 
performed and the results programmed into SSPR. An equation was derived 
that results in effectiveness values with a +1^ degree of accuracy for 
sink to base temperature ratios under a prescribed boundary. 


2.0 DISCUSSION 


Steady state analysis of a radiating fin can be accomplished 
by the use of the radiating fin effectiveness where the effectiveness is 
defined by: 

_ Actual heat radiated 

Heat radiated at the fin base temperature 


where: 




(Q/A) rei 

( a Ct^^ - T 
b s 


- radiating fin effectiveness 
(Q/A)reJ = heat flux radiated by fin 


(Eq. 1) 


T^ = fin base temperature 

T = equivalent sink temperat\ire of the environment 

S ", 

€ = emissivity 


cr = Stefan-Boltzman constant 

For a heat absorbing fin such as a solar absorber the values 
of 7] are needed for T greater than 1.0. The procedure for obtaining 

these is discussed befow. 


The fundamental differential equation for the temperature 
distribution of a radiating one dimensional rectangxilar fin can be derived 
(Ref. IB ) and is given by: 


d T _ € a 

dx^ k S 


(T^ - T 

s 


(Eq. 2) 
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where: 


T . = fin temperature which is a function of x 

X = the distance along the fin perpendicular to the edge 
which is at T = T, 

U 

5 = the fin thickness 

The value of rj can he determined as a function of T /T^» 

Nc by integration of the basic differential equation where Nc is given by: 

2 „ 3 
Nc = > ^ ^ ^b 

k S 

L = the length from the fin base to the end of the fin 
k = the conductivity of the fin material 
5 = the fin thickness 

This was accomplished (Ref. IB ) for values of T/T. betweem 0 
and 0.9* An exact solution can be derived for T /T, = 1.0 (Ref . 2B ) and is 
given by: ® 


n _ tanh (2 n/Nc ) 

' aVNT . 

=V(o +< Qj^)/f o 

a = solar absorptance of the fin surface 

Q = incident solar and albedo radiant flux 
s 

Q. = incident IR Radiant flux 
1 

This eq:mtion can be written in the form: 

II = - V 0 ^ - ^ - 5 ( 0 g)^(e - Qj. I 

by defining: 


(Eq. 3) 


(Eq. k) 
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Q 

s 


S D 


9 


L 





= temperature at x = 0 


Tl = temperature at x = L 
and one integration. 

From equation 1, the fin effectiveness is: 

T, = ..(£/.Al_reJ. . . 

/ 4 4 V 

e a (T ^ ) 

D S 

Using equation 4, this equation can be reduced to: 

<T- 




Vn 7 (1 - 9 


- 5(eg)^(i 


- 9r) (Eq. 5) 


where: 


Nc = 

k S 


The values of rj were calculated as a function of Nc , and 0 by a step wise 
integration of equation (4) to determine values of Qjj (as a function of 
Nc and 0^ ) and substituting these values into equation (5)- 

If 0j^ is assumed to equal 0^ equation (5) reduces to: 

= Nc (f- 9 M V (i+el - 50s + 1) (Eq. 6) 

S 

Equation (6) was found to result in effectiveness values 
with a degree of accuracy for sink to base temperatxare ratios under a 
prescribed boundary. This boundary is shown in Figure 1-B along with 
values of fin effectiveness. 
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